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ABSTRACT 

We report on-off pointed MAMBO observations at 1.2 mm of 61 5piteer-selected star-forming 
Q . galaxies from the 5piteer Wide Area Infrared Extragalactic Legacy survey (SWIRE). The sources 

^ ' are selected on the basis of bright 24 /im fluxes (/24pm>0.4mJy) and of stellar dominated yl&zi- 

^2 ■ infrared spectral energy distributions in order to favor z~2 starburst galaxies. The average 

' 1.2 mm flux for the whole sample is 1.5±0.2mJy. Our analysis focuses on 29 sources in the 

Lockman Hole field where the average 1.2 mm flux (1.9±0.3mJy) is higher than in other fields 
(l.l±0.2mJy). The analysis of the sources multi-wavelength spectral energy distributions indi- 
. cates that they are starburst galaxies with far-infrared luminosities from 10^^ to 10^^ '^ Lq, and 

' stellar masses of ~0.2-6x lO^^M©. 

^SJ . Compared to sub-milhmeter selected galaxies (SMGs), the SWIRE-MAMBO sources are 

among those with the largest 24/im/millimeter flux ratios. The origin of such large ratios is 
investigated by comparing the average mid-infrared spectra and the stacked far-infrared spectral 
energy distributions of the SWIRE-MAMBO sources and of SMGs. The mid-infrared spectra, 
00 , available for a handful of sources, exhibit strong PAH features, and a warm dust continuum. 

' The warm dust continuum contributes to ~34% of the mid-infrared emission, and is likely asso- 

ciated with an AGN component. This constribution is consistent with what is found in SMGs. 
The large 24/im/1.2mm flux ratios are thus not due to AGN emission, but rather to enhanced 
PAH emission compared to SMGs. The analysis of the stacked far-infrared fluxes yields warmer 
' dust temperatures than typically observed in SMGs. Our selection favors warm ultra-luminous in- 

frared sources at high-z, a class of objects that is rarely found in SMG samples. Indeed SMGs are 
not common among bright 24 /.tm sources (e.g. only about 20% of SMGs have /24^im >0.4mJy). 
Our sample is the largest Spiteer-selected sample detected at millimeter wavelengths currently 
available. 

Subject headings: Galaxies: formation - Galaxies: starburst - Galaxies: high-redshift - Cosmology: 
observations - Millimeter 
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1. Introduction 

After the discovery by IRAS of local Ultra- 
Luminous InfraRed Galaxies (ULIRGs, with hpjji 
> 10^^ L(7)) and fol l ow-up studies (see rev i ews b y 



Sanders fc Mirabell Il996l : iLonsdale etld] l2006l ). 



the advent of sensitive sub-mm bolometer ar- 
rays heralded a revolution in our understanding 
of star formation at high-red shift. Extragalactic 
surveys with SCUBA/JCMT ( Hughes e t al.lll998l; 



Scott et al.ll2002tlMortier et al.li2005; Coppin et al. 
20061 ) aiid the Ma x Planck Mill i meter Bolometer 



(MAMBO) array (jKrevsa et a l."l998^ at the In 
stitut de Radioas tronomie Mil l imetri q ue (IRAM) 
30m telescope (|Greve et al. 120041: IVoss et al 



[2006') show that the number of dusty, infrared- 
luminous galaxies (high-z ULIRGs, commonly 
"Sub-miUimeter Galaxies" or SMGs) at 2<z<3 
is several orders of magnitude higher than in 
the local universe and that these distant SMGs 
contribute significantly to the universal star - 
formation history (see e.g. ISmail et al.l |200J). 
SMGs appear to be "maximal starbursts" that 
form stars at the maximum global rate until 
a significant fraction of their original large gas 
mass is converted into a massive stellar system. 
Numerous follow-up studies indicate that SMGs 
are the most likely objects to trace the major 
episodes of star formation leading to the produc- 
tion of massive (> 2L*) ellipticals (jGreve et al 



2005t iTacconi et al.ll2006l ). This picture is consis 



tent with the claim that evolved massive ellipti- 
cals have formed the bulk o f their stars in short 
(<2Gyr), intense bursts ( Dunlop et al. 1996t 
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Blakeslee et~aLlll999l) . The evolution of SMGs into 
massive elliptical galaxies is often attributed to a 
negative feedback mechanism indu ced by Active 



Galactic Nuclei (AGN) activity ([Sanders et al 
19881 ). It is speculated that as the super-massive 



black hole (SMBH) in these systems reaches a cer- 
tain mass and luminosity, the AGN will heat the 
surrounding gas, and thus halt any further star- 
forming activity and SMBH growth. This scenario 
is supported by deep X-ray observations of SMGs 
that show that so me SMGs contain fast grow- 



ing S MBHs ( Alexander et al.l l2005al: iBorvs et al 
20051) . 



In this context, identifying a large sample of 
ULIRGs at the epoch where these processes were 
taking place, i.e. z >2, is particularly important 
to understand the formation of the most massive 
galaxies, investigate the contribution of SMGs in 
the star formation history of the Universe, and 
the presence and role of AGN activity. These 
goals have up to now been stymied by the lim- 
itations in mapping speed and sensitivity of the 
current sub-mm bolometers. The total number 
of known SMGs from all blank-field surveys with 
any instrument is less than five hundred. Their 
number and luminosity function are known from 
MAMBO/SCUBA counts, but remain uncertain 
at both low and high luminosities. Moreover, in 
order to quantify the relative importance of BH 
growth, stellar mass assembly, and the negative 
feedback induced into the host galaxies by these 
processes, it is necessary to separate any starburst 
and AGN power in these systems. Such an anal- 
ysis requires a multi-wavelength coverage and a 
good characterization of their spectral energy dis- 
tributions (SEDs) and spectra. For these reasons, 
numerous efforts are being carried out to collect 
large samples of ULIRGs at z ~2 with well sam- 
pled SEDs and spectra. 

2. Spitzer observations of sub-millimeter 
galaxies 

With the advent of wide area Spitzer surveys, 
new breakthroughs in the study of SMGs are 
now possible as the statistics of the high red- 
shift ULIRG population increases and the power 
sources in SMGs are quantified. In particular, 
the Spitzer Wide Area Infrared Extragalactic Sur- 
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vey (SWIRE 0; iLonsdale et 111120031 ) . the largest 
('--^SOdeg^) of the Spitzer suite of extragalactic sur- 
veys, with its sensitivity is able to detect many 
sources found in blank- field sub-mm/mm surveys. 
The SWIRE 5ct limits are 4.2, 7.5, 46, 47/iJy in 



the 4 IRAC (jFazio et all 120041 bands at 3.6, 4.5, 
5.8, and 8.0 /xm, res pectively, and 0.21 , 18, and 108 
mJy in the 3 MI PS (iRieke et a .120041) bands at 24, 
70, and 160 (jSurace et al.ll2005[ ji. Many pre- 
viously known SMGs lie in the SWIRE fields and 
a large fraction of them are detect ed in 2 or more 



IR ba nds (see statistics below and lClements et al 
20081 ). 

In order to estimate which fraction of SMGs 
would be detected by Spitzer to the SWIRE sen- 
sitivities, we gathered from the literature a sam- 
ple of 90 SMGs with publishe d Spitzer data, or 



that lie in the S^yIRE fields jGreve et all 12004 



Ivison et al.ll2004l: lEgami et al 



20041: iFraver et al 



_2004a"b': 'Chapm an et al.ll2004 iBorvs et allbOOsT 
Pope ct aL ,2006) . This sample is not meant to 
be complete or exhaustive, but it includes a large 
number of objects that are representative of sub- 
mm selected SMGs ('classical SMGs' hereinafter). 

From an analysis of the SMGs with fi.2mm > 
2.5 mJy or Ss^onm. >7mJy (22 'bright SMGs'), we 
expect the large majority to be detected in at least 
2 SWIRE bands. Ten (45%) of these 22 'bright 
SMGs' are detected above the SWIRE 24^m --5- 
6cr flux density hmit {^250fi3y), 7 of them (32%) 
are brighter than 400/iJy, 20 (91%) are detected at 
both 3.6 and A.bfim above the SWIRE 5cr limits, 
18% are detected at 5.8^m and 27% at 8.0/in|l. 

The density of bright (/i.2mm > 2.5 mJy or 
Sssou T n >7mJy) SMGs is ^ 230 per sq. deg. (jVoss et al, 
2006t ICoppin et al1l2006h . which means that in 
~50 sq. deg. SWIRE would detect with IRAC 
about ~10,000 sources that are typical of those 
found in sub-mm surveys, and ^^5,000 at 24/im 
(>3,000 brighter than AOOfiJy), which is 10 times 
larger than all existing MAMBO-SCUBA surveys 
put together. Appreciable SWIRE detectabil- 
ity rates are also found for the sub-mm fainter 
SMGs (3.2<S850Mm<7mJy): 36% detectable at 



^http: / /swire. ipac.caltech.edu/swire/swire. html 

^The given SWIRE 5a limits are valid for the LH field. They 

can differ by up to a factor of f .2 in the other fields. 
■^The Spitzer 70 and 160/im imaging sensitivity is such that 

extremely few high redshift ULIRGs are detectable in large 

area surveys such as SWIRE. 



24/zm (14% brighter than 400/iJy); 86% at 3.6 
and A.bfiui, 5% at 5.8/xm and 9% at 8/im, so many 
thousand more sub-mm fainter SMGs are also to 
be found among SWIRE sources, including many 
24/im-detected sources. Moreover, Spitzer offers 
the opportunity to investigate possible selection 
biases inherent to sub-mm surveys in favor of ob- 
jects dominated by cool dust. 

The identification of high-z ULIRGs among the 
large SWIRE population of galaxies is not trivial, 
however, because they are faint in the optical and 
their optical-MIR SEDs can be similar to those of 
more quiescent star forming galaxies. Moreover, 
only long wavelength data (far-IR and mm/sub- 
mm) and reliable redshifts can confirm the ex- 
treme luminosities of the ULIRG candidates, but 
only 0.4% of the SWIRE population is detected 
at A > 70/im, and most of these ULIRG candi- 
dates are at z<l or too faint in the optical to be 
foUowed-up spectroscopically. 

X-r a y, radio, IR and CO studies ([Alexander et al 
|2005at [Chapman et al] l2005t iGreve et al.l 12001" 



Val iante et al.l 2007t iMenendez-Delmestre et aL. 
.2007i r indicate that the far-IR (FIR) emission of 
the bulk of SMGs is dominated by starbursts, but 
a large fraction of them also host weak AGNs. 
Spitzer holds a key to this question because the 
diagnosti c provided by the M IR broadband SEDs 
and IRS ( Houck et al.l 2004 ) spectra can discrim- 
inate warm AGN-dominated emission. Indeed, 



using IRS spectra, iPope et al.l (|2007r ) constrain 



the AGN contribution to the MIR emission of 
SMGs to <30%. However, Spitzer alone cannot 
determine the FIR luminosity or dust tempera- 
ture of these systems because most of them are 
not detected at 70 and 160/im. 

Therefore, millimeter or sub-millimeter obser- 
vations are needed to identify these objects and 
directly measure the long wavelength emission and 
estimate their FIR luminosities. The combination 
of SWIRE and MAMBO offers the possibihty to 
explore the role of ULIRGs at z ~2 and determine 
their power sources. 

We have thus undertaken a systematic observ- 
ing program with MAMBO at 1.2 mm to search for 
luminous ULIRGs among SWIRE sources. The 
results from the first set of observations from this 
program are presented here. 

This work is organized as it follows. § |3] de- 
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scribes the selection of the best candidates from 
the four northern SWIRE fields. Their observa- 
tions with MAMBO are presented in § HI The 
spectral energy distributions (SEDs) of all sources 
and photometric redshift estimates are presented 
in § [5] Millimeter and 24 fj,m fluxes, broad-band 
SEDs, and colors of the selected sources are com- 
pared with those of sub-mm selected galaxies in 
§ [6l FIR luminosities and star formation rates 
(SFRs) are estimated by fitting greybody mod- 
els with dust temperatures typically measured in 
SMGs, and presented in § [7l Since the major- 
ity of the sources are not detected with MIPS at 
70 and 160 /zm, we measure average values from 
stacked MIPS images. The stacked FIR fluxes at 
70 and 160/Ltm, combined with the average 24 /im 
and 1.2 mm fluxes, are used to characterize the 
FIR emission of these sources, and to compare 
it with the FIR emission of SMGs in § [71 Stel- 
lar masses, near-IR (NIR) luminosities and star- 
formation timescales are discussed in § [S] The 
presence of AGN activity is investigated using dif- 
ferent tracers, MIR spectra. X-ray, and radio ob- 
servations in § m The possible nature of the 
selected sources and their relation with classical 
SMGs are discussed in § [TUl Finally, our re- 
sults are summarized in § [11] Most of the analy- 
sis described here is applied to the sub-sample of 
sources from the Lockman Hole field. We only re- 
port the observations and resulting data on the 
sub-samples from the other fields. This choice 
is due to the fact that the Lockman Hole sub- 
sample is larger and better observed than those 
in the other fields. The sub-samples from the 
other fields do not allow us to carry out an accu- 
rate statistical analysis. We adopt a flat cosmol- 
ogy with Ho = 71kms~^Mpc -\ fJM = 0.27 and 



vious works t argeting extreme 24 um sources (e.g 



Houck et al.1 l2005t lYan et all l2005 h. For 



= 0.73 (|Spergel et al.ll2003l ) 



3. Sample selection 

The most sensitive of the Spitzer MIPS bands 
for identification of distant ULIRGs is the 24/im 
band. Therefore, we have searched for extreme 
ULIRGs among SWIRE sources with /24 > 
250/zJy, corresponding to >5cr. Most systems 
above this /24 limit and z>l are expected to be 
ULIRGs. In addition distant ULIRGs will be 
faint at optical wavelengths, therefore we limited 
our overall z>l ULIRG samples to sources with 
r'>23. This is a brighter cut than adopted in pre- 



purpose, we prefer to adopt a brighter optical 
limit because a fainter cut can reject ULIRGs 
with colors like some local well known exam- 
ples, e.g.; an Arp220-like system (z=0.018 and 
L(IR) = 1.45xlO^^Lo) at z=l with an observed 
/24 = 250/iJy would have a luminosity of Lm = 
lO^^'^L© and an observed magnitude r'=23.4, 
while an Arp 220-like system with Lm = IO^'^Lq 
at z=l would have an observed r'=22.4. 

The selection based on bright 24 ^m and faint 
r'-band fluxes yields consequently sources with 
large f24fj.m/fr ratios, typically >100. Bright 
24/^m Spitzer samples with large IR/optical flux 
ratios (i.e., /24^m//i? >500) te nd to be doin i nated 



by AGN light i n the MIR faouck et al.' [20051; 



Yan et al.[ [20051: JSaiina et all ^007; Bra nd et al 



20061 : iDev et al.l I2008I ). Moreover MAMBO 
observations of these have shown that they 
are weak at 1 . 2 mm compared to sub-mm se- 



lected SMGs ([Lutz et al.[ [20051) . An additional 
selection criterion was thus applied to disfa- 
vor AGN-dominated and favor star-forming- 
dominated systems. The additional selection 
criterion is based on the evidence in the IRAC 
bands for the ~ 1.6/im rest-frame Planck spec- 
trum peak from evolved stars. Note that IRAC 
was in part designed for photometric selec- 
tion of galaxies at z=1.5-3 disp l aying this fea - 
ture ([Simpson &: Eisenhardt[ll999l: [Sawicki[[2002l ). 
Such a peak is not observed in templates of lo- 
cal AGN-dominat ed sources such as Mrk231, 
IRAS 192 54-7245 ([Berta et al.[ [2003l). or in type 



1 QSOs (jElvis et al.[[l994[ ). On the other hand, 
the SED of sources displaying such a peak can 
be well fitted by local star formation -dominated 



temp lates, such as M82, and Arp 220 ([Silva et al 



19m . 



Furthermore, we have selected sources where 
this peak lies in the 5.8/xm band in order to fa- 
vor sources at 2=1.5-3. We refer to such sources 
as '5.8 /im-peaker' (or 'b3' for bump in the third 
IRAC channel). Formally a source is defined 
'5.8 /xm-peaker' if satisfies the following condition: 
/a. 6 < /4.5 < fb.s > fs.o where / is the fiux den- 
sity in /ijy and /s.o can also be an upper limit 
(typically at ba). 

We have previously demonstrated that this se- 
lection method is effective in finding starburst- 
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dominated objects at z =1.5-3 (IWeedman et al 
20061 : iBerta et al.ll2007bHFarrah et al.ll2008l ). IRs" 



MIR spectroscopic observations of a similar sam- 
ple of 32 SWIRE sources displaying a rest-frame 
NIR bump peaking at 4.5 ^m dem onstrate that 



they are star-forming dominated (jFarrah et al 
20081 ). Most of their spectra exhibit promi- 
nent PAH emission features consistent with a 
starb urst origin of their M IR emission. Simi- 
larly, l^edmanetiD dlooi) find that all SWIRE 
sources with a NIR bump at 5.8 /xm, display PAH 
emission features in their IRS spectra. Moreover, 
8 of the sources selected here were observed in 
the MIR by the IRS i nstrument on Spitzer (6 in 



the LH; Weedman et al. 2006, and 2 in the EN2 



field, Lonsdale, in prep.), and their spectra display 
strong PAH features as expected from starburst- 
dominated emission. The IRS observations yield 
spectroscopic redshifts that, in all but one case, 
confirm the '5.8 /xm-peakers' photometric redshift 
selection method. The agreement between pho- 
tometric and IRS redshift s is di scussed further in 

2008h . Based on these 



[5] (see also iFarrah et al 



considerations and studies, we infer that the com- 
bination of a strong 24 /xm flux and a peak in the 
IRAC bands is a good indication of the presence 
of PAH emission at z ~2, and thus of a significant, 
and likely dominant, starburst contribution to the 
NIR and MIR emission. Although in some cases 
an AGN component in the MIR can be present, 
its contribution is never dominant compare d to 
the starburst component (jBerta et al.ll2007al ). 

In order to further increase the chance of select- 
ing z ~2 star-forming galaxies, we derived pho- 
tometric redshifts for all source candidates (for a 
detailed description of the procedure see § O . 

We have selected in this manner 1023 '5.8 /im- 
peaker' ULIRG candidates with /24 > 250/iJy, 
r'>23, and photometric redshift in the range 1.5- 
3, from the SWIRE's largest field (lldeg^), the 
Lockman Hole (LH). We restricted the search to 
the region in the LH with available optical imag- 
ing (~9.5deg^). At similar redshifts and for fixed 
SED shape the most luminous objects will be the 
brightest in flux. Therefore, to cull from this large 
sample the most luminous candidates we selected 
the brightest 24 /im sources, i.e. /24 >400^Jy, 
among all '5.8 /im-pcakers', and with >5ct de- 
tections in the flrst 3 IRAC bands. There are 
520 sources with /24 >400/iJy, corresponding to 



55 sources per square degree. We have also re- 
quired unambiguous (unique or not detected coun- 
terpart) source identification at all optical-MIPS 
wavelengths in the images, so that source confu- 
sion would not affect the SED interpretation. Our 
selection is likely to strongly favor sources with 
strong 7.7/im PAH features in 24/im band in the 
z=2 range. 

Finally we made estimates of the expected 
1.2 mm flux densities by fitting templates of var- 
ious local starbursts and ULIRGs to the optical 
(Ugriz) and infrared (3.6-24/im) bands. Since 
most of the candidates are not detected at 70 
or 160/im, their FIR through mm spectral shape 
is unconstrained. The MAMBO flux predictions 
therefore rest on the adoption of a long wavelength 
SED template. To select appropriate templates 
we build the rest-frame SEDs of -30 SMGs with 
known redshifts and available IRAC and 24/im 
fluxes, and normalized them at rest-frame l.Gfim. 
We then selected templates of local starburst- 
dominated ULIRGs which pass approximately 
through the midst of the distribution of 850/im 
and 1.2mm points for these objects, as shown in 
Figure [Ij Most of the SMGs have FIR SEDs simi- 
lar to Arp220 or IRAS 22491-1808 (122491 here- 
inafter). The main difference between the Arp220 
and 122491 templates, a starburst dominated and 
a composite starburst -I- AGN system, respectively, 
resides at optical and MIR wavelengths. However, 
since with some extinction the optical emission of 
122491 resembles that of Arp 220, and in the MIR 
we do not have enough spectral coverage to favor 
one template with respect to the other, we de- 
cided to fit each candidate source with an Arp 220 
template. In addition we decided to also fit a 
template with a much lower FIR/NIR flux ratio 
than Arp 220, to obtain minimum estimates of the 
predicted 1.2mm flux for each source. For this we 
chose NGC 6090, a local starburst with a similar 
SED shape to the well known starburst galaxy 
M82. Examples of these two flts are shown for 3 
sources in Figure [2l where it can be seen that the 
two templates produce similar photometric red- 
shifts from the IRAC data, but the predictions for 
the bolometric luminosity, and for the 1.2 mm flux 
density differ by a factor of about 10. 

To select the best candidates for observations 
with MAMBO, we averaged the predicted 1.2mm 
flux density for each source from both templates 
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Fig. 1. — Rest-frame SEDs normalized at 1. 6 iiia of sub- mm selected galaxi es from the l iterat ure (open 
circles: HDFN sources from IPope et al. ( 2006[ ). crosses: GOO DS sources from Borvs at al. ( 2005h . squares: 
Lockman Hole sourc es, and tr iangles: ELAIS-N2 sources from Chaoman et al. ' ("2004!), diamonds: Lockman 
Hole sources from I Greve et al.l (12004 ) . stars: FLS sources from .Fraver et al.. (,2004a3 . and plus signs: GOODS 
sources from lEgami et al.l "( 20041 ). The curves represent various starburst galaxy templates: M82 (sohd red), 
NGC6090 (dotted green), Arp220 (dashed blue), and 122491 (dot-dashed magenta). 
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densities, Arp 220 (solid curve), NGC 6090 (dotted curve), and 122491 (dashed curve). Downward arrows 
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and ordered the sources according to decreasing 
predicted averaged 1.2 mm flux. While this proce- 
dure may not predict the 1.2mm flux of an indi- 
vidual source accurately, because the actual long 
wavelength SED is unconstrained, it does provide 
a robust way to prioritize them by decreasing the 
expected 1.2 mm flux under the assumption that 
they have similar SED shapes, and taking into ac- 
count the sources redshift. 

The same selection procedure was applied to 
other three SWIRE fields, XMM-LSS (XMM), 
ELAIS-Nl (ENl) & ELAIS-N2 (EN2). A descrip- 
tion of the fields size s and multi- wavelen gth cover- 
age can be found in Surace et al.l (|2005l ). The top 
61 sources from the 4 SWIRE fields, (LH, XMM, 
ENl & EN2), were observed with MAMBO at 
1.2 mm. The largest number of sources was se- 
lected in the LH field, 29 sources, compared to 17 
in ENl, 9 in EN2, and 6 in XMM. The predicted 
1.2 mm fiuxes (averaged on all templates) range 
from 1.4 to 25 mJy. The list of selected sources, 
coordinates and optical-infrared-mm fluxes are re- 
ported in Tables [1] and m 

Note that because of the uncertainties on the 
IRAC fluxes, especially at 5.8 and 8.0 /im, which 
are less sensitive bands than those at shorter wave- 
lengths, a '5.8 /zm-peaker' {z ~2.2-3.2) can be- 
come '4.5 /Ltm' or '8.0 /im-peakers' ('b2', or 'b4') in 
different versions of the data processing. Indeed, 
some of the sources selected here, that were pre- 
viously classified '5.8/j,m-peakers' based on a pre- 



vious version of the SWIRE catalog used during 
the selection process, show a peak at shorter ('b2'; 
z ~1.5-2.2) or longer wavelengths ('b4'; z >3.2) 
or a fiat IRAC SED ('fl') in the latter SWIRE 
catalogs that is used in this work. This change is 
observed in 14 out of 61 selected objects (see last 
column of Table [1]) . 

4. MAMBO observations and results 

Observations were spread over the pool ob- 
serving sessions at the IRAM 30m telescope in 
FaU/Winter 2005/2006, using the 117 element ver - 
sion of the MAMBO array (jKrevsa et all Il998l ) 
operating at a wavelength of 1.2mm (250 GHz). 
On the first run in late October - early November 
2005 we observed most of the Lockman sample 
(29 sources), and the 6 XMM sources. During the 
Winter run, mainly early March 2006, the observa- 
tions of the Lockman sample were completed and 
exploratory observations of the 17 ENl and 9 EN2 
sources were carried out. 

We used the standard on-off photometry ob- 
serving mode, chopping between the target and 
sky at 2 Hz, and nodding the telescope every 10 
or 20 s. On-off observations were typically ob- 
tained in blocks of 6 scans of 16 or 20 lOs-subscans 
each, and repeated in later observing nights. The 
atmospheric transmission was intermediate with 
T(1.2mm) at zenith between 0.1 and 0.4. The ab- 
solute flux calibration was established by obser- 
vations of Mars and Uranus, resulting in a flux 
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calibration uncertainty of about 20%. 

On average, the noise of the channel used 
for point-source observations was about 35-40 
mJy/v^/beam, where t is the exposure time in 
seconds, consistent with the MAMBO time esti- 
mator for winter conditions. In order to identify 
a maximum number of strong sources with flux 
density ~3-5mJy within the allocated time (57h, 
i.e. ~38h of integration time), we adopted a two- 
step procedure: a first exploratory phase generally 
aiming at an rms ~1.0-1.2mJy; then longer inte- 
grations on sources showing hints of detection in 
order to reach rms « 0.6-0.7 mJy and achieve 3- 
5(7 detections of sources stronger than ^2-3 mJy 
However, we were able to apply such a complete 
procedure only to the LH sample (29 sources), 
which had the best selection and was the only one 
observable in good conditions (during night or 
early morning) both in Fall and Winter runs. The 
observations of the sources of the three other fields 
(ENl and EN2, XMM; 32 sources) were more or 
less limited to the exploratory phase. 

A major change occurred in the 30m-telescope 
control system in 2005. MAMBO pool observa- 
tions in our October-early November session were 
carried out with a well known system which had 
been in use for years. A new system, just in- 
stalled before the Winter pool session, brought 
some problems such as stronger accelerations of 
the telescope during nodding with possible ef- 
fects on tracking accuracy and temperature stabil- 
ity of the helium cryostat in some circumstances; 
anomalous noise in some MAMBO channels for a 
few cases; etc. These problems induced some de- 
lays and loss of time, but they were practically 
all well under control for most of our observations 
in February-March 2006. The data were reduced 
with standard procedures in the MOPSIC pack- 
age developed by R. Zylka. Great care was taken 
to check the quality of the data obtained with the 
new control system, and we discarded the data 
from about ten scans displaying obvious problems 
such as anomalous noise in some channels, as well 
as ~3% of the data taken in periods when we are 
not sure that the problems with the new system 
were well under control. 

Table [T] lists the resulting 1.2 mm flux densities 
and their statistical uncertainties. In the more 
completely studied Lockman sample (29 sources. 
Tabled]), more than 40% of the sources (12) show a 



signal at >2.5cr, and fluxes >2mJy (correspond- 



ing to ~5 mJy at 850 /im for z ~2, see iGreve et al 



2004). Eight of them have 1.2mm fluxes in the 
3-5 mJy range, with S/N>4-5, which guaran- 
tees the robustness of the detections. However, 
it is possible that the flux densities of the de- 
tected sources are slightly overestimated by the 
Malmquist bias, by an amount that probably does 
not exceed 10% on average. The average flux den- 
sity of the 29 LH sources is 1.88±0.28mJy (see 
Table [3]), corresponding to ~5mJy at 850 /im at 
z ~2). 

As quoted, the remaining 32 sources, from ENl 
& EN2 and XMM, were not observed as com- 
pletely. However, the average strength of the 
sources in the other fields is significantly weaker 
(1.14±0.20mJy) than in the Lockman sample. 
The fraction of sources showing at least a 2.5(t sig- 
nal is only 25% (8 sources), there is practically no 
source stronger than 3 mJy, and the average of the 
23 other undetected sources is only 0.60±0.17 mJy. 
Possible reasons for such a behavior are either a 
pointing problem caused by the system changes 
implemented in 2005, or a bias introduced by the 
slight different selection criteria (see ? 16. ip . 

Table [3] summarizes the information about the 
number of sources showing a signal at various lev- 
els and the average values for two sub-samples and 
the totality of the 61 sources. 

5. Photometric redshifts and SED fits 

In order to classify the spectral energy distri- 
butions (SEDs) and estimate photometric red- 
shifts of the sources in the sample, optical and 
IR data (<24/im) are combined and fitted with 
various galaxy templates . The SEDs are fitte d 
using the Hyper-z code ( Bolzonella et al. 2000() . 



Hyper-z finds the best-fit by minimizing the 
derived from the comparison of the observed 
SED and expected SEDs at various redshifts, 
and using a template library and the same pho- 
tometric system. The effects of dust extinction 
are taken into account by reddening the refer- 
ence templates according to a selected reddening 
law. We use the prescription for e xtinction mea- 
sured in high-redshift starbursts by ICalzetti et al. 
(I2OOOI) . We limit the additional extinction Ay 



to be less than 2.0 mag and include templates 
of highly extinguished objects. We used two tem- 
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plate libr aries, one with 105 tem plates of starburst 
galaxies ( Chary fc Elbad 2001 ). and one with 8 
templates including 2 late-spirals, 5 starbursts, 
and 1 composite (starburst+AGN) template cov- 
ering the wavelength range between lOOOA and 
1000/im. The starburst templates correspond 
to the SEDs of NGC6090, N GC6240, Arp220 
IRAS 20551-4250, and 122491 (jSilva et al.lli99i 
llBcrta 2005). The composite (AGN-^SB) template 
corresponds to the SED of the Seyf ert 2 galaxy 
IRAS 19254-7245 South (119254; iBerta et al ' 



20031 ). A description of the library a nd of the 



method can be found in lPolletta et al.l (|2007f ). 

We fit the SEDs using the data up to 24/im, 
including upper limits. Although the optical up- 
per limits place some constraints on the acceptable 
fits, for these optically faint objects the photo- 
metric redshifts are based primarily on the IRAC 
bands. The photometric redshifts are therefore 
limited in accuracy by (1) the small number of 
available photometric bands (3-4 IRAC bands, 
MIPS [24], and up to a maximum of 5 optical 
bands among Ugriz), and (2) the fact that the 
main spectral feature available for fitting (the peak 
in the stellar spectrum) is broad and detected 
through broad filters. Therefore our photomet- 
ric redshifts are expected to have uncertainties 
of ±0.5 in redshift or higher, even for the high- 
est SNR sources. In order to illustrate the range 
of acceptable solutions and photometric redshifts, 
we report multiple solutions if they give similar 
X^- The fits were also done without including the 
24/im measurements. Since the comparison with 
spectroscopic redshifts was slightly better when 
those data were included, we decided to adopt 
the results obtained using also the 24/im data. 
The inclusion of 24 /tm data in the fitting proce- 
dure can yield better results, if templates includ- 
ing dust emission, like the libraries use d here, are 
employed (see e.g. IPoUetta et al. 2007 ). 

The reliability and accuracy of the photomet- 
ric redshifts can be estimated using the available 
spectroscopic data. Spectroscopic redshifts from 
the IRS ar e available for 8 sourc es and are listed 
in Table l4l ( Weedman et al. 20061 . Lonsdale et al., 
in prep). We define the fractional error Az, the 
systematic mean error Az, the la dispersion cr^, 
and the rate of catastrophic outliers, defined as the 
fraction of sources with lAzI > 0.2. Az is defined 



as: 



and 



Az 



^phot 



1 + 



Zphot 



spec 



'^spec 



/N 



(1) 



(2) 



1 ~^ ^spec 

with N being the number of sources with spectro- 
scopic redshifts. The systematic mean error, Az, 
is 0.03. The rms of Az, (7^ is 0.08, and there 
are no outliers with |Az| > 0.2. The comparison 
between the spectroscopic and final photometric 
redshifts is shown in Figure [3] for the 8 sources. 
There is a good agreement for 5 sources, and for 
3 sources there is a difference >0.05 in Log(l-l-z). 
In order to better asses the reliability of our esti- 
mates, we applied the same technique to a sample 
of 24 fim detected SMGs from the literature with 
known spectroscopic redshift and with at least 3 
detections in the IRAC bands consistent with the 
presence of a peak at either 4.5 /im or 5.8 /im. We 
found 15 sources that satisfy these criteria out of 
the 90 SMGs for which we collected data from the 
literature (see ii l6.ip . Note that most of the sources 
from the literature do not benefit from observa- 
tions from as many bands as the SWIRE-MAMBO 
sample. After combining the SWIRE-MAMBO 
and the literature sample, Az=0.03, az=0-19, and 
we find 2 outliers (|Az| > 0.2), corresponding to 
~9% of the sample. The results obtained for the 
literature sample confirm the accuracy of our es- 
timates on the SWIRE-MAMBO sample. 

SED fits for all the SWIRE sources in the sam- 
ple are provided in the Appendix and the best-fit 
templates and photometric redshifts are listed in 
Table 12 

The redshift distribution of the sample, based 
on the photometric redshifts (solid line) and on 
the spectroscopic redshifts (dotted line) when 
available, is shown in Figure [H The distri- 
bution peaks at z ~2.0, which corresponds 
to about the minimum redshift for a 5.8 /im 
peaker. The z peak of our sources is slightly 
lower than the redshift distribution of sub-mm- 
select e d sources ( F'raver et al.l 2004bt Egami et al 



20041: iBorvs et a,l.l 120051: IChapman e t all bOOsF 



Greve et al.ll2004l) . Note that the difference could 



be even larger because of the possible overesti- 
mate of Zpfiot in case, for example, of a shift of 
the peak of the stellar bump due to an AGN con- 
tribution (e.g. a 'b2' can appear as a 'b3' and 
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Fig. 3. — Comparison between photometric and spectroscopic redshifts from IRS (Lonsdale et al., in prep. 
Weedman et al. l2006l) of the 8 SWIRE/MAMBO sources with available spectroscopic redshifts (diamond and 
triangles). Dotted lines represent 20% agreement in (I+2). Full symbols represent sources with /i.2mm>2o', 
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other three SWIRE fields. Crosses represent 15 sources from the literature sample of SMGs with available 
spectroscopic redshifts, a 24 /im detection, and exhibiting a peak in the IRAC bands. 
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its redshift be overestiinated by Az =l, see e.g. 
Berta et al.ll2007bl : lOaddi et al.H2007h . 



6. General properties of the 
MAMBO sources 



SWIRE- 



6.1. Millimeter and 24 //m fluxes 

The distribution of observed 1.2 mm flux den- 
sities (separated in sources with /i.2mm^2(T and 
<2(t) is illustrated in Figure [SJ The fraction of 
>2a sources, 41% (24% at 4tT), and the large 
average 1.2 mm flux observed in our samples 
(1.5±0.2 mjy), and especially in the more thor- 
oughly observed LH sample (1.9±0.3mJy), show 
that our criteria are quite successful in selecting 
sub-millimeter-bright galaxies from MIR samples. 

The XMM, ENl & EN2 sample show system- 
atically lower fluxes compared to the LH sample. 
Although the LH sample was more thoroughly ob- 
served than the other samples, with typical expo- 
sure times >1.5 hours, compared to ^--^40 min for 
the other sources, the origin of the different detec- 
tion rate and average fluxes is more likely due to 
intrinsic properties of the sources or of the flelds, 
e.g. cosmic variance, than to the observational pa- 
rameters. The most plausible explanation is that 
our selection introduced different biases in the se- 
lected samples possibly due to the different optical 
coverage or to the analysis carried out on the SEDs 
and on the images during the source selection pro- 
cess. The latter was indeed more extensive in the 
LH field. 

In Figure El we also show the distribution of 
1.2 mm fluxes of the 90 SMGs selected from the 
literature (see § [21) . In case of lack of flux measure- 
ment at 1.2 mm for the SMGs from the literature, 
we estimate it from the flux at 850 /xm by applying 
the following equation: 

Log{h,2mm) = -0.36 + 0.98 x Log{fs5o,.m) (3) 

where fi.2mm and /sso/im are the flux densities at 
1200 and 850 /im, respectively. This relationship 
is derived from a sample of 21 SMGs for which 
meas urements at both 850 and 1200 /xm are avail- 



able (jGreve et al.ll2004l ) 
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Fig. 5. — 1.2mm flux distribution of sources 
with /i.2mm>2cr (histogram). The mean and 
mean standard deviation of the sources with 
/i.2mm<2cr are shown as a vertical solid line and 
a gray shaded area. From top to bottom, all 
MAMBO/SWIRE sources, those in the Lockman 
Hole, those in the the other fields, and those in 
the lit erature (Fraver et al. 2004b : Egami et al 



Greve et al.l 12004 jPope et all l2006f) . 



20041: iBorvs et al.l I2OO5I: IChapman et all l20o'l 



The to- 
tal number of plotted sources is annotated 
and in parenthesis the number of sources with 

/l.2mm<2(T is givCU. 



Based on Figure 7 in Voss et al. ( 200^ and the 
850 M m source counts from SHADES ( Coppin et al.l 
20061 ). the expected number of SMGs with f 1.2mm >5 mJy, 
4mJy, and 3mJy in the 9deg^ LH field is ^-^250, 
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Fig. 4. — Redshift distribution for all observed sources. The distribution of spectroscopic redshifts is shown 
with a dotted line. From top to bottom and from left to right, all MAMBO /SWIR E sources, those in 



the Lockman Ho l e, ELAIS-Nl, ELAIS-N2, and XMM-LSS and from the l iterature (jEgami et all 12004 : 



Borvs et all l2005t IChapman etai] 120041 : iGreve et all 12004 IPope et al.|[2006f) . The number of sources is 



annotated on the upper left corner of each panel. 
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800, and 2,000, respectively. On the other hand, 
we have identified only 2, 3, and 4 sources brighter 
than 5mJy, 4mJy, and 3mJy, respectively - out 
of 29 24 bright sources in the 9 deg^ LH field. 
The small number of mm-bri ght sources in our 
sample and in the sample in iLutz et al. (l2005l) 



suggests that very luminous millimeter sources 
are not common among sources with bright 24 ^m 
fluxes (/24Mm >400^Jy). 

In Figure [HI we show the distribution of 24 ^m 
fluxes of the SWIRE-MAMBO sample and of the 
SMGs from the literature. All SWIRE-MAMBO 
sources have 24 /im greater than 400 fj,3y as im- 
posed by the selection criteria. The 24/im fluxes 
of our sample range from 0.4 to 1.5 mJy, with 
a median value of 0.76mJy and a mean value 
of 1.09±1.06mJy The SWIRE-MAMBO sources 
with /i.2mm<2CT show a similar range of 24/im 
fluxes, but fewer of the faintest 24/im sources are 
detected by MAMBO. Figure [H] clearly illustrates 
that the majority of SMGs are characterized by 
24/im fluxes smaller than 400/iJy. More specif- 
ically, out of the 90 SMGS from the literature, 
20 (22%) have 24/im fluxes greater than 400 /iJy. 
Therefore, about 78% of classical SMGs would be 
missed by our selection criteria. Indeed, more 
than half of SMGs at z ^2 are not detected at 
the SWIRE sensitivity at 24/im, --250mJy (and 
even at 5.8 or 8.0 /tm, see § ). The lowest 24/im 
flux in the SMGs reaches 21 /iJy, but there are 
16 sources that are not detected, 4 in GOODS 
(/24/jm <26-50/iJy) and the others in the SWIRE 
flelds {f24iim <250/iJy). Since our sources are 
at similar redshifts and show similar mm fluxes 
than SMGs, they either represent a subset of 24/im 
bright SMGs (~20%) or a different class of SMGs 
with enhanced MIR emission. This question will 
be further investigated in the following sections. 

6.2. Comparison with other Spitzer sam- 
ples 

A number of works have recently been pub- 
lished about the MIR properties of SMGs. They 
are based on hundreds of sources detected both 
with Spitzer from 3.6 to 24/im and with SCUBA 
at 850 /im or with MAMBO at 1.2mm (see §[2]). 
Most of these sources have been identifled in deep 
flelds mapp ed both with SCUBA or MAMBO 
and Spitzer j'Eeami et al.l[20 04: Iv ison et al.ll2004t 
Greve et al]|2004HChapman et al.ll2005HBorvs eTal 
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Fig. 6. — 24/im flux distribution of sources 
with /i.2mm>2cr (solid line) and with /i.2mm<2(T 
(dotted line). From top to bottom, all 
MAMBO/SWIRE sources, those in the Lockman 
Hole, those in the the oth er fields, and t hose in 
the l it erature ( [ Fraver et~a l. 2004b; Ega mi et al 



2004 iBorvs et al 
Greve et all 120041; 



Pope et al.l 120061) . The dashed 



20051: IChapman et al] T2004 



histogram represents the literature sources that 
are not detected at 24 /im for which we assumed 
a 5(7 upper limit (250 /iJy for 12 sources in the 
SWIRE fields, and 26/iJy for 1 GOODS SMC, 
and 506 /tJy for 3 more GOODS SMGs). The to- 
tal number of sources is annotated on the upper 
left corner in each panel. The value in parenthesis 
represents the number of sources with /i.2mm<2cr. 
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20031: IPope et all l2005l l200d llvison et all 12007 



Yun et al. 2008h . or in pointed Spitzer obser- 



vations at the position of known S CUB A or 



MAMB O sources ( Charmandaris et al .ll2004l : lFraver et al 



l2004a!B; 'Ashbv et al.' '20061; 'Valiante et al.' '200' 
Men cndcz-Dclnicstrc ct al. 2007; Pope ot al. 200 
All these studies focus on sub-mm or mm se- 
lected galaxies, while our Spitzer selection aims 
at sources bright at 24 /xm. However, it is inter- 
esting to note that our sources satisfy the selection 
criteria based on IRAC colors developed to select 
SMG ca ndidates among the IRAC population (see 
eq. 1 in lYun et al1l2008l ). 

Only a few other studies have a similar ap- 
proach to ours in characterizing the millimeter 
properties of Sp itzer mid-IR selected sources. 
Lutz et al.l (2005j study a sample of bright (/24pm ^1 mJy) 



sources with high 24/im/optical flux ratios, and 
f24:p,m/ fafim Aux ratios to favor star-forming galax- 
ies, but whi ch ended up to be l argely c ontaminated 
by A GNs (lYan et alJ l2005l l2007t ISaiina et al 



20071 ). Others study AGN-dominated s amples 



in the SWIRE survey (jPolletta et al.ll2008L Perez- 
Fournon et al., in prep.). The SWIRE-MAMBO 
sample differs from those studies because spe- 
cial care was applied to favor star forming galax- 
ies and remove AGNs that usually dominate the 
populati on of 24 ^tm brigh t and optically faint 



sources ( Brand et al. 2006f l. Since the presence 



of a bump in the rest-frame NIR is due to stel- 
lar emission, our selection disfavors sources with 
dominant AGN emission in the MIR. Our selec- 
tion thus makes this sample unique with respect 
to other 5piteer-selected samples that have been 
foUowed-up at sub-mm and mm wavelengths 

6.3. Spectral energy distributions 

The spectral energy distributions (SEDs) of our 
sources are compared with those of sub-mm se- 
lected SMGs in Figure [T] Only sources with spec- 
troscopic redshifts are shown in the right panel. 
All SEDs are shown in the rest-frame and nor- 
malized at 1.6/im in the rest frame. The SWIRE 
sources are predominantly high in the 24/im band 
at the location of the rest-frame 7.7/im PAH emis- 
sion feature, and low in the rest-frame FIR com- 
pared to the SMGs from the literature. The FIR 
emission of the SWIRE sources is in between that 
expected for the typical starburst galaxies as M 82, 
and NGC6090, and a ULIRG like Arp220, after 



normalizing them at the 1.6/im in the rest- frame. 

Figure [7] shows clearly the selection bias that 
is present in our sample: by selecting the 24/im- 
brightest sources with a peak at 5.8 /im, we have 
strongly favored sources at z ^2 with large NIR 
and MIR luminosities. The classical SMG sample 
is naturally more biased towards the most mm- 
bright systems. The location of our sources in 
Figure [7| with respect to various template tracks 
confirms that their SEDs arc intermediate between 
M 82 and Arp 220. 

6.4. Colors and nature of the sources 

We compare the 1.2 mm/24 fim flux ratio of the 
SWIRE and literature SMG samples in Figure [8] 
The SWIRE sources show on average lower flux 
ratios than the literature sources at parity of red- 
shift. The difference in the 1.2mm/24/im can not 
be simply explained by a redshift difference but is 
likely due to an excess in the MIR (i.e. 24/im) 
with respect to the FIR in the SWIRE sample 
compared to the literature sample. Such an ex- 
cess, also visible in the SEDs (see Fig. [7]), might 
be due to an AGN contribution in the MIR, or to 
enhanced PAH emission, or to a deflciency at mm 
wavelengths due to lower SFRs or higher temper- 
atures. Indeed objects characterized by warmer 
dust temperat ures are expected to show lower 
flux r atios (e.g. IChapman et al. 20031 : iLewis etliT 
2005h . 



In Figure [51 we also show for comparison 
the 1.2mm/24/im flux ratios measu red in the 
sampl e of composite sources from 'Lu tz et al 
(|2005l ). and in high-z radio galaxi es (HzRG) de- 



tecte d at both 850 /im and 24/im (iSevmour et al 
2OO7I ). The 1.2 mm flux is derived from the 



850 /j,m by applying e q uatio n [3l These two sam- 
ples, f rom iLutz et al.l ( 2005h . and Sevmour et al 



()2007l) . show even smaller 1.2 mm over 24 fim flux 
ratios. In this case, the most plausible explana- 
tion for such small ratios is the contribution from 
AGN-heated hot dust emission to the MIR since 
both samples contain mainly AGNs. Since our 
sources show intermediate flux ratios in between 
those seen in the literature SMGs and HzRGs, it 
is quite possible that an AGN contributes to their 
MIR, more than typically observed in SMGs and 
less than observed in HzRGs. 

To better investigate the origin of such an ex- 
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Fig. 7. — Rest-frame SEDs of SMGs from the literature (black symbols as in figure [T]) and of the 
SWIRE/MAMBO sources with 1.2 mm fluxes >2a (full circles, red: Lockman Hole, green: XMM-LSS, 
ELAIS-Nl, or ELAIS-N2) normalized at 1.6 /im. The 70 and 160 /xm fluxes correspond to the normalized 
stacked fluxes of all SWIRE/MAMBO sources with 1.2 mm fluxes >2a. The curves represent starburst tem- 
plates: M82 (solid red), NGC6090 (dotted green), Arp220 (dashed blue), and 122491 (dot-dashed magenta). 
All sources are shown on the left panel and only sources with a spectroscopic redshifts are shown on the 
right panel. 



cess we consider in more details the SMGs that 
fall in the same region covered by the SWIRE 
sources, and analyze the MIR spectra when avail- 
able. From the literature sub-sample with avail- 
able spectroscopic redshifts and available 24 /im 
flux (65 sources), we find 3 sources that fall in the 
same region in the /i.2mm//24/jm vs z diagram as 
our sources, J123635. 59-1-621424.1 (/i.2mm//24/im 
= 1.3 5) at z = 2.005 from the GOODS North 



and their SEDs are very different from those of our 
sources as a dominating AGN component is clearly 
seen in the MIR. 

The best and probably unique way to constrain 
the presence and the contribution of an AGN com- 
ponent in the MIR is to measure the equivalent 
width of the PAH features. A warm dust contin- 
uum due to AGN-heated hot dust will indeed de- 



field (jBorvs et al 

(/l.2mm/ /24/im 



20051) . CSMMJ163650+4057 
3.07) at z = 2.384 from the 



crease their equivalent width (see e.g.lClavel et al 



2000: Valiante et al.ll2007l: Menendez-Delmestre et al 



ELAIS-N2 field (|Greve et al.l2004[ ). and SMMJ10523 8 
(f 1.2m m / f 24am = 3.34) at z = 1.852 (jchapman et al 



20051). The first two sources are AGNs, their MIR 
SEDs are consistent with a power- law model, their 
optical spectra show broad emission lines, and one 
source is also a luminous X-ray source. The lat- 
ter source is instead classified as starburst galaxy 
based on its optical spectrum, and it also shows 
the highest 1.2 mm/24 /im flux ratio of the three 
selected sources. This analysis thus reinforces the 
hypothesis that an AGN component might be re- 
sponsible for the lower flux ratios observed in the 
SWIRE sample, but this is based on only 2 sources 



i2007t ISaiina et al.ll2007t iDesai et al.ll2007l ). This 
19+i^?iii5Jelnent can be performed using MIR spec- 
tra. 

Based on the IRS spectra of similar samples 
and the few available here, a signiflcant AGN 
contribution is disfavored for the majority of 
our sample. MIR spectra from IRS are avail- 
able onl y for 8 sources in ou r SWIRE-MAMBO 
sample (jWeedman et al. I l2006l Lonsdale et al., in 
prep.). Altogether, 6 out of these 8 sources show 
PAH dominated spectra and weak or absent con- 
tinuum, implying large equivalent widths and neg- 
ligible or absent AGN contribution. They are, in 
approximative order of PAH strength, LH-07, LH- 
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Fig. 8. — Observed flux ratios, f i. 2 mtn / f 24 fim, as a function of redshift for SMGs from the literature (black 
crosses), for SWIRE /MAM BO sources with .fi.2 mm>2a (circles: LH, triangles: other fields), and for high- 2; 
radio galaxies (blue squares ISevmour et aljEoOTt ) . Expected values for various starburst and AGN templates 
are also shown: Arp220 (sohd blue), M82 (dot-dashed purple), NGC6090 (dotted red), 119254 (dashed ma- 
genta), and QSOl (dot-dot-dot-dashed green). Filled symbols represent SWIRE/MAMBO sources with 
spectroscopic redshifts, open symbols represent sources with photometric redshifts. Redshift uncertainties 
are not reported for clarity, but they are of the order of 0.5 for the photometric redshifts, and <0.1 for 
the spectroscopic redshits. f 1,2mm fluxes for the literature sample ar e derived from /a sn ii.m. when not avail- 
able. The average 1.2mm/24/im flux ratio of 24/xm-selected galaxies ( Lutz et al.ll2005f) is shown as a dotted 
horizontal black line. 
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06, LH-26, EN2-04, LH-12, and LH-05. The only 
exceptions are 2 sources (LH-23/B6, and EN2- 
03) that show, in addition to PAH emission, a 
warm dust continuum that might be associated 
with an AGN component. We do not find any 
trend between the presence of continuum, the 
strength of the PAHs, the 1.2 mm flux, and the 
fi.2mm/ f24fj.m flux ratio. Therefore, for this anal- 
ysis we consider the average spectrum of all 8 
sources, independently of their properties. The 
average IRS spectrum of all 8 SWIRE-MAMBO 
sources is shown in Figure [H We also show for 
comparison the composite spectru m of 12 SMGs 



from the HDFN (iPope et al.ll2007). the IRS spec- 



trum , of a starburst, NGC3079 (jWeedman et al 
20051 ). a power-law model representing emission 
from warm dust, and the sum of the starburst 
spectrum and of the warm dust continuum. The 
average spectrum of the SWIRE-MAMBO sources 
clearly shows PAH features, implying a strong 
starburst component, but also a warm dust contin- 
uum. The average spectrum is well reproduced by 
the sum of the starburst spectrum and the warm 
dust continuum. The continuum is modeled with 
a power-law model (/i, oc A"^^), with spectral in- 
dex aiR=2. We have apphed the same analysis on 
the average spectrum after removing the 2 sources 
with larger continuum contamination and for the 
sources with >2a and <2a 1.2 mm fluxes and the 
result is still valid with only small differences that 
would require a larger sample to better quantify. 
The contribution of the power-law component to 
the average spectrum is ~34% in the 6-12 /xm rest- 
frame wavelength range. Assuming that the 24 /im 
flux includes a 34% contribution from an AGN on 
average, the AGN-subtracted fi.2fim/ f24tj.m flux 
ratio would be 1.5 x higher. Even after such a cor- 
rection, our sources remain systematically on the 
lower bound region of the flux ratio distribution. 
Note that for consistency we would have to ap- 
ply a similar correction to the literature SMGs as 
well. Indeed, the majority of classical SMGs also 
show both PAH features and a w eak warm dust 
continu um in their MIR spectra (IVahante et al 



2007; M enendez-Delmestre et al.ll2007tlPope et al 



2007|). The PAH equivalenth withs of our sources 
and of classical SMGs are indeed very similar, as 
shown in the left panel of Figure [31 Thus, the hy- 
pothesis that our sources might have an enhanced 
24 /im flux due to the contribution from AGN 



emission, is disfavored. However, our sources are 
more luminous at 24 /xm than SMGs. Their PAHs 
and MIR continuum are thus, on average, more 
luminous than in SMGs. 

7. Far-infrared properties 

In order to estimate FIR luminosities, and star- 
formation-rates (SFRs) of the SWIRE-MAMBO 
sources, their FIR SEDs need to be character- 
ized. Since most of the sources are too faint to 
be detected at 70 and 160/im with MIPS at the 
SWIRE sensitivity (typically 18, and 108 mJy, re- 
spectively), we can only use the observed 1.2 mm 
flux and make reasonable assumptions on their 
FIR SEDs. We assume that the FIR-mm emission 
of our sources is similar to that seen in SMGs, 
and model their FIR-mm emission with a grey- 
body model normalized at the observed 1.2 mm 
flux or 3(7 upper limit in case the flux is <2a. To 
constrain the temperature of the FIR-mm emit- 
ting dust, it would be necessary to have at least 
two flux measurements at FIR-mm wavelengths. 
Since the dust temperature can not be directly es- 
timated for our sources, we flx the temperature 
to the values typically observed in SMGs after 
taking into account the dependency of the dust 
temperature on the redshift. From the sample of 
SMGs with multiple sub- m m and mm data stud- 
ied by iKovacs et al. ( 2006h . we derive: 



Tdu.t - 3.9 + 9.5 X (1 + z) (K). 



(4) 



We us e a dust emissivity index /3=1.5 as in lKovacs et al 
(|2006l ). Higher values of ^ give slightly large r tem- 
peratures (see Fig. 3 in Kovacs et al. 20061 ). and 
larger FIR luminosities. At parity of flux and tem- 
perature, the FIR luminosity derived with /3=2 is 
on average 2.3 times larger than the FIR luminos- 
ity derived with /3=1.5. Two sub- mm/mm flux 
measurements are available only for 4 sources in 
our sample (Kovacs et al., in prep.). In these four 
cases, the estimated temperatures are consistent 
with the measured ones, giving support to our 
assumption. The dust temperatures range from 
23 K to 41 K, with a mean value of 33±4K. In or- 
der to take into account the redshift uncertainty, 
the greybody fits are performed three times: 1) 
at the redshift of the sources, Zhesti 2) at a lower 
redshift, z-=Zbest — 0.12 x (1 -f Zbest), and 3) at 
a higher redshift, z^-=Zbest + 0.12 x (1 + Zbest)- 
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Fig. 9. — Average IRS spectrum of the 8 SWIRE sources with available II IS spectra (black solid curve). 
Left panel: Also shown the composite spectrum of 12 SMGs from the HDFN (jPope et al.ll2007l l after scaling 
it to match the SWIRE average spectrum (red dashed curve). Right panel: IRS spectrum of a starburst 
(NGC 3079; green cure), a power law model with spectral index 2 representing warm dust continuum emission 
(blue curve), and the sum of the two (magenta curve). The expected locations of PAH emission features are 
labeled in both panels. 



The coefficient 0.12 corresponds to 1.5 times the 
dispersion found in the comparison between pho- 
tometric and spectroscopic redshifts (see §[5]). The 
FIR luminosities are reported in Table[5] The FIR 
luminosities for the sources with a >2a signal at 
1.2mm range from 1 to 6.3 xlO^^L©, confirming 
that our SWIRE-MAMBO sources are ULIRGs. 
This is true for the majority of the sources also at 
z — z^. A similar luminosity range is obtained for 
the rest of the sample, but they are likely upper 
limits to the true FIR luminosity. 

The greybody models are always below the 
70 /im, and 160 fim upper limits of each source, 
implying that the assumed temperatures are never 
overestimated. In case of higher temperatures, the 
FIR luminosities would be higher than estimated. 

We carry out an additional analysis to char- 
acterize the average FIR SED of our sources. We 
estimate average fluxes at 70/^m and 160/im by co- 
adding all the 70/xm and all the 160/^m images to 
create a stacked 70/Ltm image and a stacked 160^m 
image. The co- adding techniqu e used in th i s wor k 
follows the method described in iDole et al.l(|2006h . 



This technique was applied to all MIPS images 
(24, 70, and 160 fiui) and only to the sources in 
the LH field since the number of observed and de- 
tected sources is significantly larger than in the 
other fields. The stacking of the 24/im images was 
carried out to validate the method by comparing 
the stacking result with the mean flux of the single 
measurements. In the rest of the analysis, we will 
use the mean of the single measurements at 24/^m 
instead of the stacking value. Two sources, LH- 
18, and LII-28, were excluded from the stacking 
because the former is located at the edge of the 
160 /im image, and the latter is detected at 70 fim. 
We divided the remaining 27 sources in the Lock- 
man Hole sample in two sub-groups: sources with 
1.2mm fluxes above 2a (12), and below 2a (15). 

At each of the positions of the sources in each 
sub-group, we extract a 72" x 72" sub-image from 
the fuU-fleld SWIRE mosaics. The sub-images are 
placed into a stack twice, once with their orig- 
inal orientation, and again with a 90° rotation, 
to remove instrumental effects such as gradients 
affecting the stacked flux. However, since the 
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MIPS maps are alre ady filtered followin g; the tech- 
Da), the rota- 



niques described in iFraver et al. 
tion should not give a significant difference. In- 
deed, a test on the the 70 /xm stacks with and 
without images rotation yields a difference of 1.4% 
which is negligible compared to the uncertainty. 
The mean flux image is computed from the stack, 
and the flux measured in an aperture (with sky 
subtraction in an annulus included). An uncer- 
tainty estimate is derived using a bootstrap tech- 
mque (|Efron fc Tibshirani|[l99i) . in which sam- 
ples of the subimages are drawn with replacement 
from the image stack, and the mean image and 
aperture flux measurement are recomputed. The 
standard deviation of the aperture flux measure- 
ments of 200 repetitions of this process yields the 
uncertainty estimate. We have validated these un- 
certainty estimates by stacking positions randomly 
offset by 60 to 80 arcseconds from our source po- 
sitions; the standard deviation of 500 such stacks 
is slightly lower than the uncertainties from the 
bootstrap method. All of the measured stacked 
fluxes are reported in Table [H The mean of the 
single measurements at 24 /itm agrees in all three 
cases with the stacked mean 24 /im flux. 

The stack of all 27 LH sources results in a detec- 
tion at both 70 and 160/im, at ~3.8tT. The stack 
of sources >2a 1.2 mm fluxes is well {^\a) de- 
tected at both wavelengths. The stack of the sam- 
ple with i\.iram < 2cr does not result in a detection 
at either wavelength. Note that the significance 
level of these measurements reflects mainly differ- 
ences in noise level due to small number statistics. 
However, these measurements provide some con- 
straints on the average FIR flux and spectrum of 
our sources. 

The average FIR SEDs (from 24 to 1200 ^m, 
including the stacked 70 and 160 //m fluxes) of the 
LH sub-sample with /i.2mm>2o' is shown in Fig- 
ure [ini The figure also shows the predicted av- 
erage flux at 70 and 160/im assuming the grey- 
body models used to compute the FIR luminosi- 
ties. The stacked mean fluxes are well above the 
greybody models predictions (solid black curve in 
Figure fTO|l. This discrepancy cannot be explain by 
any of our uncertainty (in redshifts for example). 
At 160 /im the mean observed flux predicted by 
the greybody models is 4 times lower than mea- 
sured in the stacked images. This discrepancy sug- 
gests that our sources are more FIR luminous and 
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Fig. 10. — Average FIR SEDs of LH sources with 
f 1.2mm >2cr (black full circles). The 70 /im and 
160/im average fluxes correspond to the stacked 
values. The grey curves represent redshifted grey- 
body fits (see Eq. [3]) of all 27 LH sources with 
fi.2mm >2cr. The thick solid black line repre- 
sents the median of all the greybody fits shown 
in the Figure as grey curves. The dashed black 
line represents a greybody model fit to the average 
160 /im and 1.2mm fiuxes assuming the mean red- 
shift, 2;=2.03, and a temperature T=42.5 K. As ex- 
pected the median 1.2 mm flux of the models and 
observed are the same, but at 160/i,m, the median 
model predicts a flux that is 4.9 times smaller than 
the measured value. The triangles at 24 /im and 
850 /im represent the average observed fluxes o f 
the sample of 31 HDFN SMGs (IPope et al.ll2"006h . 
and the triangles represent the stacked fluxes at 
70 /im and 160 /im of a subsample of 26 sources 
from the same sample (jHuvnh et al.ll2007 ). 
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warmer than estimated by the greybody models 
that are based on the SMGs temperatures. In or- 
der to fit the SWIRE-MAMBO sources average 
FIR SED, the dust component should be on aver- 
age 10 K higher than predicted. Indeed, a grey- 
body model that fits the average FIR SED as- 
suming the mean redshift of the LH sub-sample, 
^=2.03, has a temperature of 42.5 K (dashed black 
curve in Figure), instead of 32.7 K as predicted by 
eq. m Note that while the excess at 160 /zm can 
be explained by a temperature effect, at 70 /im 
(>20 fim in the rest-frame) a different contribu- 
tion from small grains is more likely. In this work, 
we do not investigate the difference at 70 /xm, as 
we focus our analysis on the FIR part of the spec- 
trum. 

For comparison with the classical SMGs, we 
consider the FIR p roperties of >30 SMG s from 
the HDFN sample (|Pope et al.ll2006l l2007t ). The 
mean FIR SED of this sub-sample is also shown 
in Figure [H The 850 ^iui fluxes of these SMGs 
range from 1.7 to 20.3 mJy (corresponding to 0.7- 
8.3 mjy at 1.2 mm), and show a mean value of 
5.7mJy (corresponding to 2.4mJy at 1.2mm), 
consistent with our SWIRE-MAMBO sources with 
/i.2mm>2cr. Their average 70 and 160 /xm fluxes, 
derived fro m the stacking of a sub-sample of 26 
sources bv iHuvnh et all (^2007^ ■ are significantly 
lower than those measured in the SWIRE sample. 
Such a difference suggests that the HDFN SMGs 
are, on average, characterized by cooler dust tem- 
peratures and lower FIR lumi nosities . Note that 
the sub-sample used by Huvn h et al.l (j2Q07) does 
not include 2 sources that are detected in the 
MIPS images. However, even including those 
sources the average FIR SED of classical SMGs 
would remain well below that of our sources. 

In summary, the analysis of the average FIR 
SED of the SWIRE-MAMBO sources with /i.2mm>2cr 
indicates that their FIR luminosities are higher 
than estimated using the greybody models and 
that the FIR emitting dust is at higher temper- 
atures than assumed. On average, our sources 
are thus more luminous and warmer than classical 
SMGs. 

8. Stellar masses 

Our IRAC-based selection corresponds to a 
rest-frame NIR selection, therefore, in absence 



of a significant AGN contribution at these wave- 
lengths, we are directly sampling the stellar com- 
ponent. Moreover a NIR selection is minimally 
affected by dust extinction. Therefore it is possi- 
ble to derive accurate stellar masses, by means of 
spectrophotometric synthesis. We focus this anal- 
ysis on the better observed and larger LH sample. 
We use the c ode Sim-P hot-Spe c (SPS), devel 



oped in Padova (jBerta et al . 2004; ,Poggianti et al 
20011 ) to estimate the stellar mass in our sources. 
The code performs mixed stellar population 
(MSP) spectro-photometric synthesis. The adopted 
synthesis stellar population (SSP) library is 
based on the Pa dova evolutionary sequences 
of stellar r aodels ( Fagotto et al.l Il994a ,b) and 



isochrones ( Bertelli et al.l Il994[ ). and was com- 



puted by assuming a solar metalicity and a lSalpeter 
(jl955l ) initial mass function (IMF) between 0.15 
an d 120 M( 7) . Th e SSP spectra are built using 
the |Pickleq ( 19981 ) stellar sp ectral library, and ex- 



tended to the NIR with the lKurucd (|l993l ) stellar 
atmosphere models. Nebular features ar e added 
based on the ionization code CLOUDY ( Ferlandl 
1996). The spectra thus obtained provide a reli- 
able description of simple stellar generations up 
to ~3 /im in the rest-frame. Note that our mod- 
els do not take fully into account the behavior 
of AGB stars, in other words they underestimate 
the contribution from intermediate-age SSPs. We 
estimate that this results in overestimating stellar 
masses of at most a factor of 2 (see a discussion in 



Berta et al.ll2007a|) . Another overestimating fac- 



tor com es from the assumpti on of a Salpeter IMF , 
indeed a lKroupa et al.l (|1993I ) or a IChabrieij (|2003h 
IMF would yield about 50% lower masses (see 
iBerta et al.ll2007al) . 

The mixed stellar population synthesis com- 
bines several phases, adopting a different star for- 
mation rate (SFR) for each age. This code is well 
suited for modeling stellar populations of star- 
forming galaxies, because it allows several, in- 
dependent episodes of star formation during the 
life of a galaxy, for example resembling multi- 
ple merger-driven starburst events. Age-selective 
extinction is applied, assuming that the oldest 
stars have abandoned the dusty medium long ago. 
Since disk populations are on average affected by a 
moderate extinction (Av~<l mag, e.g. KennicuttI 
1998I ). we adopt a maximum allowed absorption 
for stars older than 1 Gyr of Ay ~ 0.3 - 1.0 mag- 
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nitudes. For younger populations the color ex- 
cess gradually increases, but is limited to Av<5. 
The best fit is sought by means of minimiza- 
tion. In case of non-detections 3tT upper limits are 
adopted. For each source in the sample, the code 
explores the SFR-extinction parameter space us- 
ing an Adaptive Simulated Annealing algorithm 
and records all the attempted models. The uncer- 
tainty in the stellar mass estimate is then given by 
all those solutions whose values lie within 3a 
from the best fi t. Details on the f itting method 
are described in lBerta et al.l (|2007ar ). 

Assuming that the total energy absorbed in 
the ultraviolet (UV)-optical domain is processed 
by dust in the thick molecular clouds embedding 
young stars and re-emitted in the MIR and FIR 
(5-1000 /im), we exploit the 24 /im observed flux 
to constrain the amount of dust and the luminos- 



ity o f young stellar populations (see iBerta et al 
2004 for more details). 



The estimated stellar masses and associated un- 
certainties are listed in Table[71 The stellar masses 
range from 2.3 to 56 xlO^°M0, with a median 
value of 1.8x10" Mq. Although the near-IR SED 
is dominated by stellar emission, it is possible that 
an AGN component is also present. Thus these es- 
timates should be considered as upper limits to the 
true stellar masses. It is also possible that the de- 
rived stellar masses overestimate the true masses 
because our models do not take into ac count the 
TP-AGB contribution (jMarastonl [20051 ) . Indeed 
such a contribution might yield lower masses, es- 
pecially at the redshifts of our sources. The stellar 
masses as a function of redshift are shown in Fig- 
ure [TT] We find that at parity of redshift, the 
stellar masses can vary by a factor of 4. No dif- 
ference in stellar mass is observed between 1.2 mm 
detected and non-detected objects. Similar masses 
are observed in objects with a wide range of lu- 
minosities and viceversa. For example, all high-z 
(z >1.5) sources with available spectroscopic red- 
shifts have similar FIR luminosities, ~ lO^^'^L©, 
but their stellar masses can vary by a factor of 5, 
i.e. M,=8-20x101°Mq. 

8.1. NIR luminosities 

In order to carry out a consistent compari- 
son between the stellar masses of the SWIRE- 
MAMBO sources with the SMGs found in the 
literature, we compare their rest-frame luminos- 



ity at 1.6 /xm, where the stellar emission peaks. 
The rest- frame luminosity at 1.6 /im is estimated 
by interpolating the observed IRAC fiuxes. The 
estimated values as a function of redshift are re- 
ported in FigurefTSl They range from 5.8x 10^° Lq 
to 3.8x10^^Lq, and have a median value of 
9.5x 10^"'^ L0. We compare the luminosities of 
sources with a spectroscopic redshift and those 
with only a photometric redshift, and sources with 
1.2 mm fiuxes greater or lower than 2a. No appar- 
ent difference in NIR luminosity is observed among 
these sub-samples. We also divide the SMGs from 
the literature in two sub-groups (shown as tri- 
angles and diamonds in Figure [T^ . those with a 
clear stellar bump in the IRAC bands and those 
for which a bump is not visible because of a dif- 
ferent spectral shape or because there are no de- 
tections at A=5.8 and 8.0 /im. There is no sig- 
nificant difference in NIR luminosities between 
the SMGs with and without a stellar bump. The 
SWIRE-MAMBO sources are systematically more 
luminous in the NIR, suggesting that they might 
be more massive than literature SMGs. On av- 
erage our sources are 3 times more luminous at 
1.6 /J,m than sub- mm selected SMGs. This result 
is not surprising since our selection is based on 
a detection at 5.8/im where the SWIRE 5a limit 
is ~50/iJy, while the literature sample is, on av- 
erage, characterized by lower 5.8 fim fluxes, and 
both samples are at similar redshifts. 

9. AGN activity 

It is well established that starburst sources 
dominate the bulk population of SMGs with 
only a minor contribution from AGN. Indeed, 
although a weak AGN emission is observed in 
a sig niflcant fraction of SM Gs, at X-ray ener- 
gies (I Alexander et al.l l2005al) . in optical spectra 



lines ('chapman et alj_ 20051). and in MIR spec- 



tra (.Valiante et a l . 2007tlMenendez-Delmestre et al 
20071: iPope et all 120071 ). its contribution to their 
bolometric luminosity is usually modest (<10- 
20%). 

In § 16.41 the contribution from an obscured 
AGN to the observed MIR emission of our sources 
was invoked to explain their unusual colors and 
SEDs. Based on the available IRS spectra, we 
estimate, on average, a 34% contribution from a 
warm dust continuum in the MIR (6-12 nui in the 
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Fig. 11. — Stellar masses as a function of photometric (open symbols) or spectroscopic (full symbols) redshifts 
for the SWIRE Lockman Hole sample (29 sources). Circles represent sources with /i.2mm>2<T and triangles 
represent sources with /i.2mrn<2cr. Uncertainties on the stellar masses correspond to 2cr. 



22 



13.0 



12.5 



O 



12.0 - 



1 1.5 - 



1 1.0 - 



10.5 - 



X SWIRE {f^,^^>2G) 
o SWIRE z/-'"" 



SWIRE z 



■phot 




O 



10.0 
0. 



A Lit. Zph„t & bump 

o Lit. Zphjjt & no bump 

^ Lit. Z3pg, & bump 

♦ Lit. Z3p3^ & no bump 
■ HzRG 



_l I I I I I I I I I I I I I I I I I I I I I I I I I I I L_ 



1.0 1.5 2.0 2.5 

Redshift 



3.0 



3.5 



Fig. 12. Monochromatic luminosity at 1.6 /im in the rcst-framc for SWIRE/MAMBO sources (circles), for 
literature SMGs (triangles and diamonds), and for high-z radio galaxies (squares). Full symbols represent 
sources with spectroscopic redshifts and open symbols with photometric redshifts. Triangles represent sources 
with a typical stellar bump in the IRAC bands and diamonds sources where the bump is not present or not 
detected. Large crosses are overplotted on SWIRE/MAMBO sources with fi.2mm>2a. 
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rest-frame). Such a component m ight be gener- 



ated by an AGN, or by H II regions ijLaurent et al 



20001 : [Popita et aLll2006l ). However, this estimate 
is based on only 8 sources. Moreover, the IRS 
targets include the brightest 24 fim sources of the 
sample and might thus be biased towards sources 
with enha nced PAH emission or larger AGN con- 
tribution ( Brand et al. hood ). To further investi- 



gate the presence of an AGN in our sample, we 
searched for possible signatures using X-ray and 
radio data. 

X-ray observations from Chandra of moderate 
depth (70 ks) are available for 3 LH sources, LH- 
10, LH-1 9, and LH-28, from t he C/iandra/SWIRE 
survey [IPolletta et a.1.1 l2006l) . and the CLASX 
survey (|Yang et all \2004 ). Although the XMM 
field has a wide shallow (10 ks) X-ray coverage 
from XMM-Newton only one sourc e, XMM-03 



has b een observed by XMM-Newton (jPierre et al 



20071 ). and it is at the edge of the field. Out 
of 4 sources with X-ray coverage, only LH-10 
is detected in the X-ray. The X-ray 0.4-8 keV 



flux of LH-10 is 2.1x10 erg s cm ^ s 



hardness ratio is HR=:— 0.7 (jYang et al 



and th e 
2004h . 



as expected for a steep unabsorbed power-law 
spectrum. The estimated 0.4-8 keV luminos- 
ity is 5.7x10'*'^ ergss""'^, consistent with being 
AGN powered and with those measured in X- 
ray de tected SMGs in HDFN ("Alex ander et al. 
2005b| ). Although, based on the hardness ratio, 



the X-ray spectrum appears unabsorbed, the de- 
rived H-iMm/^x-ray ratio, ~100, is much higher 
tha n commonly observed in AGNs, 0.3-10, (se e 
e.g. iFranceschini et alllioOSl : IPolletta et al.ll2007b . 
The large ^AdiR/hx-ray ratio suggests that the 
X-ray luminosity is underestimated and that the 
observed X-ray emission corresponds to a small 
fraction of the intrinsic X-ray emission, perhaps 
due to a reflected component (warm scattering), 
while the intrinsic X-ray emission is completely 
absorbed. The lack of X-ray detections of the 
other 3 sources does not rule out the presence 
of an absorbed X-ray luminous AGN {Lx-ray > 
lO"'^ ergss"^) or an unabsorbed AGN of moderate 
luminosity [Lx-ray < 10^^ ergss"^). Indeed the 
detection limit of the Chandra observations corre- 
sponds to X-ray luminosities >10^'* ergss"^ for an 
unabsorbed AGN at the redshift of our sources. 
An AGN of lower luminosity or absorbed, as typ- 
ically found in SMGs, would not be detected at 



our depths. 

Deep VLA radio observations at 20 cm are 
available in a 0.4 d eg^ area centered oii the Chan- 
dro/SWIRE field (|Polletta et al.ll2006l Owen et 
al., in prep.) in the LH, and shallow radio observa- 
tions are a vailable from the VL A FIRST survey in 
aU fieldfQ (jCondon et al.lll998h . LH-11, and LH- 
22 are the only sources that fall in the deep VLA 
field and are both detected, with 20 cm fluxes of 
82 and 51 /iJy, respectively. A third source, LH- 
28, is detected by the FIRST survey. The radio 
fluxes are reported in Table [TJ The radio source 
associated with LH-28 is extended, 6'.'66x3'.'16, 
corresponding to 57kpcx27kpc. The bright radio 
flux and radio extension, combined with the lack 
of mm detection, indicate that the radio emission 
of LH-28 is powered by a jet and it is, thus, an 
AGN. LH-11, and LH-22 are not resolved in the 
radio images, <1'.'24 and 0'.'87, respectively. Al- 
though their 1.4 GHz rest-frame radio luminosities 
are large, respectively, 10^*-^ and 10^''-^ W Hz~^, 
they are not in excess of what is expected based 
on the FIR-radio relationship found in both star- 
burst and radio-quiet AGNs ((Condon et al. 199^. 
Therefore, the available radio data do not imply 
or rule out the presence of an AGN in these two 
luminous radio sources. 

In summary, the available X-ray and radio data 
indicate that at least 1 out of 3 sources contains 
an AGN. This fraction is consistent with that de- 
rived from IRS spectroscopy (>30%), but all these 
samples are limited by small number statistics and 
there is no overlap among them, leaving the ques- 
tion on the AGN contribution and role in our 
SWIRE-MAMBO sample unanswered. Moreover, 
the subset of sources observed with the IRS might 
be biased towards a certain type of sources, for ex- 
ample with enhanced PAHs, and might not be rep- 
resentative of the whole SWIRE-MAMBO sample. 
Yet, since there is no apparent difference in the 
1.2 mm/24 /im flux ratios of the X-ray and radio 
AGNs and the IRS sample, the estimated AGN 
contribution might be valid for the whole SWIRE- 
MAMBO sample. However, altough it is quite 
likely that an AGN is present in a large fraction 
of our sources, it is not expected to significantly 
contribute to their MIR emission, as it is the case 



^The VLA FIRST source cata log and images are available 
from [http: / / sundog.stsci.edu| 
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for the classical SMGs (see e.g. lPope et al.ll2007l) . 
Indeed, even by taking this contribution into ac- 
count, we still find smaller /i.2/im//24/im flux ra- 
tios, larger FIR luminosities, and a FIR peak at 
shorter wavelengths than in classical SMGs. Thus, 
as already stated in § 16. 4[ it is unlikely that the 
difference between the SWIRE-MAMBO sources 
and classical SMGs is due to a more important 
AGN contribution in our sample. 

10. Discussion 

In the previous sections we have investigated 
the multi-wavelength properties of the SWIRE- 
MAMBO sources selected here and compared 
them with those of classical SMGs. Our sources 
are characterized by redshifts, 1< z <3, and 
mm fluxes in the same range as classical SMGs. 
Compared to classical SMGs, they are in general 
brighter at mid and near infrared wavelengths. 
Their FIR luminosities, constrained for the sources 
detected at 1.2 mm, are in the ULIRG range, 
>1O^^L0. Their FIR emission is likely domi- 
nated by warmer dust than estimated in classical 
SMGs, but this result needs to be confirmed by 
additional sub-mm data. Their stellar masses are 
large, 10^°-^~"-^ Mg, and their NIR luminosities 
and, likely stellar masses, are also higher by about 
a factor 3 than in classical SMGs. Although some 
of these estimates are characterized by large un- 
certainties, it is a clear that our sources do show 
some differences compared to classical SMGs, or 
that they belong to a subclass of SMGs with larger 
MIR/FIR flux ratios than the bulk of SMGs (Fig- 
ure [8]). Indeed, SMGs with a relatively small 
1.2nim/24/xm flux ratio represents only a minor- 
ity of all SMGs reported in the literature. 

The possible origin of the observed differences 
between our sources and classical SMGs are dis- 
cussed here. Warm dust emission is an indi- 
cator of strong overall heating of the dust. A 
large infrared luminosity and faint optical/UV 
emission implies elevated dust opacity, which 
in turn might imply a larger density of inter- 
stellar dust closer to the heating sources, or a 
more efficient heating mechanism relative to the 
dust (jPale et all 120071 : iKovacs et al.ll2006h . Dust 
heating is predominantly fueled by the emission 
from massive (>8Mf:^) st ars in star-forming re- 
gions ( Kovacs et al.l l2006l ) . or by an AGN radia- 



tion field. The peak of IR emission of cold dust can 
thus shift towards shorter wavelengths in case (1) 
the radiation field g ets higher because of younger 
stellar populations (jPiovan et al. l2006t) . or (2) a 
significant contribution from AGN-heated dust is 
present. Thus, the difference in the FIR emis- 
sion peak between our sample and classical SMGs 
(see § [7]) might be an indication of more massive 
stars relative to dust mass, or of a more signifi- 
cant AGN contribution in our SWIRE-MAMBO 
sources compared to classical SMGs. 

The latter scenario has been investigated by 
measuring the intensity and equivalent width of 
the PAH features. Indeed, weak PAH emission 
are indicative of an important AGN contribution. 
Since, the average MIR spectrum, and thus the 
PAH features equivalent width, of our sources and 
of classical SMGs are consistent (see Figure [9|), 
the idea of having a significant AGN contribution 
is disfavored. This result supports the alternative 
scenario of having a higher contribution of massive 
stars relative to the dust mass. However, the avail- 
able data do not allow us to test this hypothesis 
and rule out other explanations, perhaps related 
to the star fo rming regions geom etry and distribu- 



tion (see e.g. iFarrah et al. I l2008h . 



Our results indicate that our selection has re- 
vealed a class of ULIRGs that is under-represented 
in current SMG samples. The SWIRE-MAMBO 
sources indeed cover a different region in the 
luminosity-temperature {L — T^) diagram than 
classical SMGs and bridge t he gap between the lat- 
ter and local UL IRGs (see IChapman et al. 20031 
Lewis et al.ll2005h . Similar results have been ob- 
tained by a study on a similar sample of Spitzer- 
selected sources with known spectroscopic red- 
shifts from strong PAH features, and multiple 
FIR detections (Younger et al., in prep.). 

11. Summary and conclusion 

Using the SWIRE survey and the MAMBO 
camera on the IRAM 30m telescope, we have iden- 
tified 21 ULIRGs at z ~l-3 with average 1.2mm 
fiuxes of 2.9±0.2mJy. They are characterized 
by starburst-dominated SEDs and by unusually 
large MIR/FIR fiux ratios. These sources were 
discovered among a sample of 61 ULIRG z >1 
candidates selected in the SWIRE survey that 
were followed-up with MAMBO observations at 
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1.2 mm. 

The primary selection criteria of the SWIRE 
sample of ULIRG candidates are a peak in the 
IRAC 5.8/^m band due to the red-shifted near- 
infrared spectrum of evolved stars, a bright 
(>400/xJy) detection at 24/^m, and a faint op- 
tical counterpart (r'>23). The average 1.2 mm 
flux of the whole sample is 1.5±0.2mJy. 

The optical-MIR SEDs of the selected sources 
are consistent with those of starburst galaxies. 
Photometric redshifts are estimated by fitting the 
SEDs with galaxy templates. For eight sources, 
spectroscopic redshifts from MIR spectra obtained 
with IRS on board of Spitzer are available and are 
in reasonable agreement with the estimated pho- 
tometric redshifts. 

Our analysis focuses on 29 sources in the Lock- 
man Hole field where the average 1.2 mm flux 
(1.9±0.3mJy) is higher than in the XMM-LSS, 
ELAIS-Nl & ELAIS-N2 fields (l.l±0.2mJy). 
However, we also present basic results for 32 
sources from the other fields, which were less 
thoroughly observed. The fraction of LH sources 
with 1.2 mm fluxes >2a is ~40% (24% at 4cr; 
see Table [3]) . This fraction is much higher than 
previously found for Spitzer MIR selected sam- 
ples (23%, 18 %, and 10% at >2, 3, and 4cr, 



respectively; iLutz et al.l 120051 ). We attribute 



this difference to the fact that earlier samples 
favored AGN-dominated, rather than starburst- 
dominated systems. Our sample, on the other 
hand, shows systematically lower 1.2mm/24/im 
flux ratios than the majority of Spitzer-detected 
sub-millimeter-selected SMGs. Such a difference 
is attributed to a bias in our selection in favor of 
24/^m-bright systems compared to the selection of 
SMGs. 

Far-infrared luminosities, and star formation 
rates are estimated assuming a 2;-dependent dust 
temperature as observed in classical SMGs, and 
normalizing a greybody model to the observed 
mm flux. Average 70 and 160/im fluxes from 
stacked Spitzer/MIPS images are measured. The 
mean FIR stacked fluxes are ~4 times larger than 
the fluxes predicted by the greybody models, sug- 
gesting higher FIR luminosities and a FIR peak 
at shorter wavelengths than in classical SMGs. 
The SWIRE-MAMBO sources are ULIRGs in the 
2;=1.5-3 range with FIR luminosities from 10^^ 
to 10^'^'^Lq. Stellar masses, estimated by model- 



ing the optical-infrared data with stellar synthesis 
populations, are also large, '--^O.2-6xlO^^M0, sug- 
gesting that these systems might be the precursors 
of the most massive ellipticals. 

Compared to sub-millimeter selected galaxies 
(SMGs), the SWIRE-MAMBO sources are among 
those with the largest 24/im/millimeter flux ratios. 
The origin of such large ratios is investigated by 
comparing the average mid-infrared spectra and 
the stacked far-infrared spectral energy distribu- 
tions of our sources and classical SMGs. 

The mid-infrared spectra, available for a hand- 
ful of sources, show strong PAH features and a 
warm dust continuum. The warm dust continuum 
contributes to ^34% of the mid-infrared emission, 
and is likely associated with an AGN component. 
This constribution is consistent with what is found 
in SMGs, and indeed the PAH equivalent widths 
are, on average, similar in our sources and in 
SMGs. Thus, even taking into account the AGN 
component, the 24y^m/1.2mm flux ratios remain 
larger than typically observed in SMGs. The hy- 
pothesis of contribution from AGN emission to the 
MIR is thus disfavored, and it is more likely that 
the large ratios are due to more luminous PAH 
emission that in classical SMGs. 

The analysis of the stacked far-infrared SEDs 
yields warmer dust temperatures than typically 
observed in classical SMGs. Thus, our selection 
has favored a warm class of ultra-luminous in- 
frared sources at high-z that is rarely found in 
current SMG samples. Indeed classical SMGs 
are not common among bright 24 /im sources as 
the SWIRE-MAMBO sample (e.g. only about 
20% of SMGs have f2i^m >0.4mJy). Our sam- 
ple occupies a region in the L — Td diagram that 
is under-populated by classical SMGs, implying 
that ULIRGs at high-z might exhibit a similar 
range in dust tempera t ures as loca l ULIR Gs (see 



Chapman etal\ l2003t iLewis et al.l l2005h . The 



analysis of our sample suggests that warmer 
ULIRGs sre brighter MIR sources and might ex- 
hibit more luminous PAH emission than cooler 
ULIRGs. However, to confirm such a correlation, 
an analysis of the MIR spectra for a larger sample 
of 24 /xm-selected objects with known mm fluxes 
would be necessary. This analysis will be soon pos- 
sible with other Spitzer samples that have been or 
will be observed with the IRS and with MAMBO 
(Younger et al., in prep., Lagache et al., in prep.). 
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Our sample is the largest Spitzer-sclcctcd sam- 
ple detected at millimeter wavelengths currently 
available. 
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A. SED fits and photometric redshifts 

The best and acceptable fits to the optical (Ugriz) and MIR (3.6-24 /im) SEDs of each source, obtained 
as described in §[5] are shown in Figures [T51 [Til [T51 and \W[ We report a maximum of three solutions, those 
with the lowest 3 x^. In case a spectroscopic redshift is available, the best-fit template at the spectroscopic 
redshift is reported as dotted curve. Although the templates extend to FIR wavelengths, the FIR data are 
not used in the fitting procedure and in the computation. The parameters (template, and redshift) of 
each solution are reported in Table 2] 

Eight sources have spectroscopic redshifts and are shown in Figures [T3] and [151 It can be appreciated from 
this figure that with only the broad NIR bump to constrain the photometric redshift, using the broad IRAC 
filters, a redshift accuracy of 10% in {l+z) is a reasonable expectation. In some cases, a NIR bump can be 
present in sources at lower redshift. Indeed the presence of warm dust continuum can make th e NIR bump 
or the turnover at Sum less pronounced and even shift the stellar peak to longer wavelengths (Bert a et al 



2007bl : Daddi et al. 2007 ). Also strong emission from the 3.3/im PAH feature falling into the 5.8 /im filter can 



produce a peak in the IRAC SED that would be wrongly interpreted if no optical or NIR data are available. 
These two effects increase the uncertainty associated with the photometric redshift estimate. Some examples 
can be seen in the fit to source LII-12 (zspec=l-89 and Zphot—'^-Q^) and LII-23 (zgpec^O.QS and Zphot=^-21) 
in Figure [H 
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Fig. 13. — Optical-IR SEDs of MAMBO sources in the Lockman Hole field. The solid curves represent 
different fits to the optical-MIR data (up to 24 /im). Blue, red and green curves corrrespond to fits with 
increasing x^. The template name and photometric redshifts of each fit are annotated. The dotted curve 
corresponds to the best-fit template plotted at the spectroscopic redshifts, The source ID number, and the 
spectroscopic redshift when available, are reported on the upper left corner in each panel. Downward arrows 
represent 5a upper limits at optical and infrared wavelengths, and 3ct upper limits at 1.2 mm. 
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Fig. 13. — Continued 
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Fig. 13. — Continued 
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Fig. 14.— Optical-IR SEDs of MAMBO sources in the ELAIS-Nl field. Symbols as in Figure [HI 
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Fig. 14. — Continued 
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Fig. 15.— Optical-IR SEDs of MAMBO sources in the ELAIS-N2 field. Symbols as in Figure [HI 
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Fig. 16.— Optical-IR SEDs of MAMBO sources in the XMM-LSS field. Symbols as in Figure [131 
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Table 1 

Basic properties of the SWIRE-MAMBO sample 



Source ID 


lAU Name* 


0:2000 52000 


z 


f 1.2mm 


fro 


/160 /20cm 


Notes'' 






(degrees) 




(mJy) 




(mJy) 





Lockman-Hole 

5(7 



LH- 


■01 


SJ105007.26+571651.0 


162.53026 


57.280842 


1.35 


5.65+0.74 


<18 


<108 




b3 


LH- 


■02 


SJ103639.57-f575346.6 
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5.29+0.85 


<18 


<108 




b3 


LH- 
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b3 


LH- 
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4(7 
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LH- 
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161.50743 


59.154690 


1.66 


2.06+0.79 


<18 


<108 


0.082+0.004 


b3 


LH- 


■12 


SJ104839.33+592149.0 


162.16386 


59.363609 


1.89^= 


2.56+1.08 


<18 


<108 




b4/IRS(B8) 


LH- 


■13 


SJ104656.24+594008.0 


161.73434 


59.668900 


Icr 

1.70 


1.92+0.98 


<18 


<108 




b2 


LH- 


■14 


SJ104724.25+572150.6 


161.85106 


57.364052 


2.03 


1.95+1.03 


<18 


<108 




b3 


LH- 


■15 


SJ104822.05+583828.2 


162.09186 


58.641171 


2.80 


1.13+0.77 


<18 


<108 




b4 


LH- 


■16 


SJ103809.63+591650.1 


159.54013 


59.280590 


2.40 


1.57+1.08 


<18 


<108 




b4 


LH- 


■17 


SJ105750.75+571652.7 


164.46147 


57.281300 


1.85 


1.47+1.10 


<18 


<108 




b3 


LH- 


■18 


SJ105943. 13+580849.1 


164.92970 


58.146969 


2.16 


1.46+1.09 


<18 


<108 




b3 


LH- 


■19 


SJ104716.02+593201.3 


161.81674 


59.533691 


2.11 


1.87+1.58 


<18 


<108 




b2/X 


LH- 


■20 


SJ104816. 14+592943.3 


162.06725 


59.495350 


1.90 


1.04+1.00 


<18 


<108 




b2 


LH- 


■21 


SJ103704.32+584755.1 


159.26801 


58.798641 


<lcr 
2.63 


0.75+0.90 


<18 


<108 




b3 


LH- 


■22 


SJ104646.02+590517.5 


161.69173 


59.088188 


1.54 


1.40+1.78 


<18 


<108 


0.051+0.004 


b2 


LH- 


■23 


SJ103809.17+583226.1 


159.53822 


58.540581 


0.98^= 


0.93+1.17 


<18 


<108 




b3/IRS(B6) 


LH- 


■24 


SJ103647.44+591819.5 


159.19765 


59.305408 


1.25 


1.19+1.57 


<18 


<108 




b3 


LH- 


■25 


SJ104059.86+574536.7 


160.24942 


57.760201 


2.13 


0.74+1.06 


<18 


<108 




b3 


LH- 


■26 


SJ103856.94+585244.0 


159.73723 


58.878891 


1.88^= 


0.44+0.93 


<18 


<108 




b3/IRS(B7) 


LH- 


■27 


SJ104553.76+565930.9 


161.47398 


56.991909 


2.01 


-0.10+1.45 


<18 


<108 




b3 


LH- 


■28 


SJ104942.11+561234.5 


162.42545 


56.209579 


2.82 


-0.14+1.31 


19+1 


<108 


2.73+0.14 


b4/X 


LH- 


■29 


SJ103500.98+573847.5 


158.75407 


57.646519 


2.16 


-1.49+1.18 


<18 


<108 




b3 



ELAIS-Nl 



Table 1 — Continued 



Source ID lAU Name' 



ENl 


01 


SJ160343 


08+551735 


8 


ENl 


02 


S J 160440 


45+543103 





ENl 


03 


SJ160651 


22+545338 





ENl 


04 


SJ161658 


14+535319 


3 


ENl 


05 


SJ161412 


14+541927 


2 


ENl 


06 


SJ160502 


00+540639 





ENl 


07 


SJ160409 


72+541808 


9 


ENl 


08 


S J 162002 


45+542933 


8 


ENl 


09 


SJ161950 


53+553740 


9 


ENl 


10 


SJ160945 


58+534757 


9 


ENl 


11 


SJ161632 


58+540736 


9 


ENl 


12 


SJ161358 


92+542812 


2 


ENl 


13 


SJ155904 


74+551137 


2 


ENl 


14 


SJ161753 


53+541840 


9 


ENl 


15 


SJ161642 


78+542355 


6 


ENl 


16 


SJ161414 


48+554255 


4 


ENl 


17 


SJ160611 


83+551832 


5 


EN2 


01 


S J 163734 


44+415151 


4 


EN2 


02 


SJ163528 


58+403843 


9 


EN2 


03 


S J 164058 


63+415432 


2 


EN2 


04 


S J 163022 


17+404957 


5 


EN2 


05 


SJ163506 


01+411735 


1 


EN2 


06 


SJ163748 


16+404922 


1 


EN2 


07 


SJ163729 


80+415722 


1 


EN2 


08 


S J 163424 


64+410954 


5 


EN2 


09 


S J 163755 


89+413416 


8 



XMM-01 SJ021926.24-045212.9 
XMM-02 SJ022521.66-053545.7 



02000 ^2000 2 f 1.2mm. 

(degrees) (™Jy) 



/to /i60 /20cm Notes'" 
(mJy) 



240.92949 
241.16856 
241.71341 

244.24226 
243.55058 

241.25833 
241.04051 
245.01019 



55.293282 


2.37 


2.81+0.78 


<17 


<104 


b3 


54.517490 


1.93 


2.46+0.81 


<17 


<104 


b3 


54.893879 


1.90 


2.11+0.69 


<17 


<104 


b3 


2(T 












53.888691 


2.00 


2.82+0.99 


<17 


81+4 


b3 


54.324211 


1.88 


2.09+1.03 


<17 


<104 


b3 


la 












54.110840 


2.37 


1.15+0.86 


<17 


<104 


b3 


54.302460 


1.40 


1.25+1.01 


<17 


<104 


b3 


54.492710 


2.03 


1.26+1.02 


<17 


<104 


b3 



<lo- 



244.96054 


55 


628040 


1 


69 


1.07+1.73 


<17 


<104 


b2 


242.43991 


53 


799412 


2 


38 


0.74+1.23 


<17 


<104 


b4 


244.13573 


54 


126919 


1 


55 


0.56+1.04 


<17 


<104 


b3 


243.49551 


54 


470051 


2 


11 


0.04+1.81 


<17 


55+4 


b3 


239.76977 


55 


193661 


1 


95 


-0.11+1.37 


<17 


<104 


b3 


244.47304 


54 


311352 


2 


39 


-0.27+2.00 


<17 


<104 


b3 


244.17825 


54 


398769 


1 


96 


-0.46+1.99 


<17 


<104 


b3 


243.56035 


55 


715401 


2 


22 


-0.99+2.26 


<17 


<104 


b3 


241.54930 


55 


309040 


2 


15 


-0.61+1.38 


<17 


<104 


b3 






ELAIS-N2 

















30- 












249.393494 


41 


864269 
la 


1 


85 


2.52+0.70 


<18 


<124 


248.869080 


40 


645519 


1 


57 


2.60+1.47 


11+1 


110+4 


250.244308 


41 


908958 


1 


77c 


1.15+0.80 


<18 


<124 


247.592380 


40 


832649 

<1(T 


1 


69"= 


1.17+0.83 


<18 


<124 


248.775055 


41 


293072 


2 


06 


1.10+1.16 


<18 


<124 


249.450653 


40 


822811 


1 


95 


0.90+1.21 


<18 


<124 


249.374176 


41 


956150 


2 


15 


0.35+1.55 


<18 


<124 


248.602676 


41 


165131 


1 


96 


0.08+1.07 


<18 


<124 


249.482895 


41 


571320 


2 


18 


-0.24+2.02 


<18 


<124 






XMM-LSS 














4cr 












34.85934 


-4 


870250 

3cr 


1 


89 


3.20+0.71 


<24 


<126 


36.34024 


-5 


596040 


1 


93 


2.45+0.61 


<24 


103+4 



b3 
b2 

fl/IRS 
b4/IRS 

b3 
b3 

b3 
b3 



b3 
b3 



Table 1 — Continued 



Source ID 


lAU Name'' 


O2000 
(dej 


fcooo 
frees) 


z 


f 1.2mm 

(mjy) 


/to 


/l60 /20cm 
(mJy) 


Notes'^ 


XMM-03 


SJ021933.90-043335.5 


34.89124 


-4.559850 
Icr 

-3.388180 

<lcr 

-6.004930 
-5.522100 


2.11 


2.17±0.71 


<24 


<126 


b3/X 


XMM-04 


SJ021925.60-032317.4 


34.85665 


2.12 


1.46±0.80 


<24 


<126 


b3 


XMM-05 
XMM-06 


SJ021947.11-060017.7 
SJ022512.20-053119.6 


34.94630 
36.30082 


1.01 
1.01 


0.78±0.90 
0.90±1.10 


<24 
<24 


<126 
<126 


b3 
b2 



Note. — The sources are ordered by field and signal-to-noise at 1.2 mm. 

^'SJ stands for SWIRE J. SWIRE JHHMMSS.ss+DDMMSS.s is the official lAU source name for sources discovered 
in the SWIRE fields. 

''Shape of the IRAC SED (b2: peak at 4.5 /im, b3: peak at 5.8 /^m, b4: peak at 8.0 Mm. fi: flat IRAC SED). Sources 
with IRS observations are marked with IRS. B# refers to the source ID in lWeedman et al. ! ||2006|) . Sources with X-ray 
coverage are marked with X. 

'^Spectroscopic z from IRS spectrum. IRS spectra of LH sources are published in lWeedman et"ar 1 ll2006t) . those of 
EN2 wiU be published in a future publication (Lonsdale et al., in prep.). 

"^X-ray source 425 in CLASX catalog l|Yang et al.ll2004^ . F(0.4-8 keV) = 2.1x10-^^ erg cm-^ s~i and HR = -0.7. 



Table 2 

Optical & infrared data of the SWIRE-MAMBO sample 



Source ID 


u 


g' 


r' 
(Vega) 




z 


/3.6 


/4.5 


/s.s 
(/*Jy) 


/s.o 


/24 


T w ni 
J_in-U-L 




ZO.'il 


9*^ 7"^ 


99 07 






y4 


1 OS 


^A7 
<^4 { 


( OO 


T W 09 








"~~-»9'< 

^zo.o 




41 


6\0 
OZ 


79 
1 Z 


^A7 
<.4 ( 


oO / 


T w m 




9/1 fi9 


9/1 C^/1 


9Q fll 
Zo. Oi 




oy 


C^9 
OZ 


7ft 
1 o 


04 


791 
( Zl 


ijri-U'± 






OA C\f^ 


\9Q ^ 




/I /I 


DO 


1 on 

lUU 


QQ 

yy 


707 


T w nc; 

ijil-UO 




9/1 AC\ 


9Q /I 

ZO. 4l:0 


99 7Q 
ZZ. / O 




77 


yo 


QQ 

yo 


7a 


1 i=;o9 


T w r\fi 




9"=; t^o 


--i.9A 4 

^Z'rt.'Tl 


9^ fll 
Zo. Ol 




Do 


Q9 

yz 


Q7 

y I 


79 
1 z 


1 079 


T PT 07 




9"=; 9*^ 


OA fiS 


9*^ QS 

ZO. yo 




•^7 
O ( 


Al 
4 ( 


fiO 
DZ 


AO 
4Z 


1 90*^ 
IZUO 


Lrl-Uo 






-^9/1 /I 


\9Q 




DO 


fi/1 
o4 


liD 


Di 


f^l 
OiO 


T HQ 




94 Sfi 


9zL 9*^ 


9^? 




(^0 
DU 


1 D 


yo 


oO 


ooZ 


ijrl-lU 




9(^ "^Q 


9f^ U.A 


OA 11 




9S 
Zo 


40 


DO 


^ Al 
^4 ( 


4o4 






91^ fli^ 


9*^ 70 
Zo. 1 U 


0'\ 1 7 
Zo. i I 


ZZ.OD 


70 
( U 


Q1 


yo 


70 


1 Vi^A 
lUo4 








99 (^9 
ZZ.OZ 






1 


yo 


1 90 
iZU 


1 9c; 
IZO 


1 /1 1 O 

i4iy 


J_jrl-lo 


0/1 on 


9/1 fiQ 
Z4.00 


--^9/1 A 


>Zo.O 




f^O 


fi7 


^1 
/ 


^ Al 
<4 / 


/ 00 


LH-14 




25.33 


24.36 


23.24 




42 


53 


84 


<47 


795 


LH-15 


>24.3 


>25.2 


23.97 


23.33 




54 


75 


114 


123 


638 


LH-16 




>25.2 


>24.4 






36 


46 


62 


81 


1318 


LH-17 






24.88 


23.81 




81 


102 


104 


68 


982 


LH-18 




>25.2 


>24.4 


>23.5 




20 


28 


59 


<47 


967 


LH-19 


>24.3 


>25.2 


>24.4 


>23.5 




33 


40 


57 


<47 


488 


LH-20 


>24.3 


>25.2 


>24.4 


>23.5 




35 


44 


48 


<47 


474 


LH-21 




>25.2 


23.58 


>23.5 




54 


77 


103 


104 


733 


LH-22 


>24.3 


>25.2 


>24.4 


>23.5 


>23.6 


31 


40 


38 


36 


446 


LH-23 




>25.2 


>24.4 


>23.5 




26 


33 


40 


<47 


1111 


LH-24 




24.77 


24.07 


23.00 




75 


88 


103 


81 


899 


LH-25 




>25.2 


>24.4 


>23.5 




31 


42 


51 


46 


430 


LH-26 




>25.2 


>24.4 


>23.5 




26 


45 


47 


<47 


1105 


LH-27 




>25.2 


23.51 


>23.5 




43 


54 


64 


<47 


673 


LH-28 






23.65 






30 


40 


46 


53 


499 


LH-29 




>25.2 


>24.4 


>23.5 




31 


39 


58 


<47 


451 



ENl-01 


>23.4 


>24.94 


>24.04 


>23 


18 


>21 


9 


35 


48 


63 


59 


698 


ENl-02 


>23.4 


>24.94 


>24.04 


>23 


18 


>21 


9 


37 


53 


68 


44 


771 


ENl-03 


>23.4 


>24.94 


>24.04 


>23 


18 


>21 


9 


41 


52 


57 


46 


823 


ENl-04 


>23.4 


24.42 


>24.04 


>23 


18 


>21 


9 


39 


49 


77 


37 


657 


ENl-05 


>23.4 


>24.94 


>24.04 


23 


14 


>21 


9 


30 


40 


53 


<47 


607 


ENl-06 


>23.4 


>24.94 


>24.04 


>23 


18 


>21 


9 


27 


37 


<40 


53 


458 


ENl-07 


>23.4 


>24.94 


>24.04 


>23 


18 


>21 


9 


54 


67 


62 


<47 


560 


ENl-08 


>23.4 


>24.94 


>24.04 


>23 


18 


>21 


9 


19 


32 


44 


34 


628 


ENl-09 


>23.4 


>24.94 


>24.04 


>23 


18 


>21 


9 


70 


93 


88 


56 


581 


ENl-10 


>23.4 


>24.94 


>24.04 


>23 


18 


>21 


9 


24 


37 


56 


51 


446 


ENl-11 


>23.4 


>24.94 


>24.04 


>23 


18 


>21 


9 


37 


<8 


51 


67 


673 


ENl-12 


>23.4 


>24.94 


>24.04 


>23 


18 


>21 


9 


44 


53 


80 


63 


879 


ENl-13 


>23.4 


>24.94 


>24.04 


>23 


18 


>21 


9 


35 


44 


51 


46 


694 


ENl-14 


>23.4 


>24.94 


>24.04 


>23 


18 


>21 


9 


35 


46 


75 


72 


550 


EN 1-1 5 


>23.4 


>24.94 


>24.04 


>23 


18 


>21 


9 


33 


42 


55 


<47 


631 


ENl-16 


>23.4 


>24.94 


>24.04 


>23 


18 


>21 


9 


28 


35 


48 


<47 


546 


ENl-17 


>23.4 


>24.94 


>24.04 


>23 


18 


>21 


9 


26 


41 


65 


43 


599 
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Table 2 Conl iiuKd 



Source ID 


u 


9' 


r' 


i' 


z 


/a. 6 


/4.5 


fs.s 


fa.o 


/24 








(Vega) 










(MJy) 






EN2-01 


>23.4 


>24.94 


>24.04 


>23.18 


>21.9 


45 


58 


59 


40 


814 


F,N2-n2 

i_Jl ^ ^ \J ^ 


>23.4 


>24.94 


>24.04 


>23.18 


>21.9 


36 


53 


48 


58 


928 


EN2-03 


>23.4 


>24.94 


>24.04 


>23.18 


>21.9 


49 


54 


67 


53 


1237 


EN2-04 


>23.4 


23.70 


>24.04 


>23.18 


>21.9 


83 


115 


129 


173 


1521 


EN2-05 


>23.4 


>24.94 


>24.04 


>23.18 


>21.9 


37 


51 


68 


37 


765 


EN2-06 


>23.4 


>24.94 


>24.04 


>23.18 


>21.9 


29 


48 


60 


43 


733 


EN2-07 


>23.4 


>24.94 


>24.04 


>23.18 


>21.9 


29 


38 


56 


<45 


472 


EN2-08 


>23.4 


>24.94 


>24.04 


>23.18 


>21.9 


34 


44 


57 


42 


685 


EN2-09 


>23.4 


>24.94 


>24.04 


>23.18 


>21.9 


19 


31 


44 


35 


421 


XMM-01 


>25.0 


>24.0 


>24.2 


>24.2 


>23.6 


47 


72 


99 


74 


1277 


XMM-02 


26.08 


23.86 


24.45 


23.84 


22.98 


36 


47 


58 


<66 


1512 


XMM-03 


>25.0 


24.39 


24.36 


22.32 


21.41 


70 


77 


93 


58 


1194 


XMM-04 






>24.2 






47 


63 


113 


95 


1253 


XMM-05 


26.12 


25.70 


24.86 


24.20 


22.98 


54 


62 


59 


50 


1398 


XMM-06 


>25.0 


>24.0 


>24.2 


>24.2 


>23.6 


62 


71 


<58 


<66 


1109 



Note. — Typical uncertainties are ~5% for the optical and the IRAC data, and 10% for MIPS. 
Upper limits correspond to 90% completeness in the optical and to 5cr in the infrared. 



Table 3 

Statistical information on MAMBO L2mm results 



Sample 


N 


>2cr 


N 

>3cr 


>4mJy 


< h.2rara{AU) > 

(mjy) 


< /l.2mm(> 2(7) 

(mJy) 


< /l.2mm(< 2(t) 

(mJy) 


Total 


61 


21 (34%) 


16 


3 


1.49±0.18 


2.90±0.22 


0.75±0.14 


LH 


29 


12 (41%) 


9 


3 


1.88±0.28 


3.19±0.36 


0.95±0.22 


EN1+EN2+XMM 


32 


9 (28%) 


7 





1.14±0.20 


2.51±0.12 


0.60±0.17 



Note. — Uncertainties correspond to the standard deviation of the mean. 
''Mean 1.2 mm flux density of sources with f 1.2 mm > 2<t. 
''Mean 1.2mm flux density of sources with f 1.2mm < 2a. 
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Table 4 

Best-fit templates and photometric redshifts 



oource iU 


z 


A 

Ay 


lemplate 


LH-01 


1.345 


1.43 


SBScy 


LH-02 


1.925 


1.69 


cc_96 




1.934 


1.95 


N6090 


LH-03 


2.674 


1.04 


122491 




2.332 


0.26 


cc_81 


LH-04 


2.943 


1.04 


120551 


LH-05 


1.583 


0.39 


ce_86 




1.734 


1.04 


spil 4 


LH-06 


1.791 


0.78 


ce_76 




1.803 


1.56 


122491 


LH-07 


2.225 


1.17 


122491 


LH-08 


2.014 


0.78 


CG_66 




2.035 


1.04 


Arp220 


LH-09 


2.283 


0.91 


120551 


LH-10 


2.234 


1.43 


120551 


LH-11 


1.661 


0.52 


cc_105 




1.848 


1.43 


spi2c_3 


LH-12 


1.783 


0.00 


cc_81 




2.652 


0.65 


122491 


LH-13 


1.699 


0.26 


cc_96 


LH-14 


2.032 


1.17 


122491 




2.158 


0.39 


cc_81 


LH-15 


2.804 


0.39 


cc_81 




3.235 


0.78 


120551 


LH-16 


2.406 


0.13 


119254 




2.450 


1.95 


122491 


LH-17 


1.850 


0.52 


N6090 


LH-18 


2.164 


1.95 


N6090 


LH-19 


2.111 


1.82 


N6090 




2.359 


1.04 


120551 


LH-19 


2.362 


0.65 


ce_86 


LH-20 


1.896 


1.43 


NfiOQO 




1.990 


0.52 


ce.76 




1.945 


0.52 


ce_96 


LH-21 


2.630 


0.39 


cc_81 




2.915 


0.65 


120551 


LH-22 


1.535 


1.04 


cc_76 




1.831 


0.78 


122491 


LH-23 


1.207 


1.95 


cc_86 


LH-24 


1.254 


1.43 


SBSey 




1.756 


1.17 


122491 


LH-25 


2.133 


1.43 


120551 




2.180 


1.56 


N6090 


LH-26 


1.853 


1.95 


ce.81 


LH-27 


2.014 


0.52 


122491 




2.108 


0.00 


ce.81 


LH-28 


2.815 


0.65 


122491 


LH-29 


2.164 


1.69 


N6090 




2.342 


1.04 


120551 
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Table 4 — Continued 



„ in 
oourcs lU 






J-GHipicltS 


ENl-Ol 


2.369 


1.69 


120551 


ENl-02 


1.934 


1.95 


N6090 




1.999 


1.43 


ce_91 


ENl-03 


1.896 


1.95 


N6090 




2.002 


0.91 


ce.86 


ENl-04 


2.002 


0.91 


122491 




2.136 


0.13 


ce_81 


ENl-05 


1.879 


0.26 


ce.81 


ENl-06 


2.366 


1.95 


N6090 




2.699 


1.17 


120551 


ENl-07 


1.397 


1.95 


ce.71 




1.836 


0.65 


N6090 


ENl-08 


2.026 


1.95 


N6090 




2.127 


1.95 


ce.81 


ENl-09 


1.685 


0.91 


N6090 




1.867 


0.78 


cc_66 


ENl-10 


2.379 


1.95 


N6090 




2.744 


1.30 


120551 


ENl-11 


1.545 


1.95 


120551 




3.500 


0.00 


N6240 


ENl-12 


2.111 


1.95 


N6090 




2.299 


1.30 


ce_81 


ENl-13 


1.945 


1.95 


N6090 




2.164 


1.04 


ce.81 


ENl-14 


2.392 


1.95 


N6090 


ENl-15 


1.960 


1.95 


N6090 




2.205 


0.91 


ce.81 


ENl-16 


2.221 


1.04 


ce.81 




2.594 


1.04 


122491 


ENl-17 


2.152 


1.95 


N6090 




2.342 


1.95 


120551 



EN2- 


■01 


1, 


.848 


0, 


.78 


ce.91 






1, 


.969 


1, 


.04 


122491 


EN2- 


■02 


1, 


.570 


1, 


.95 


ce.96 


EN2- 


■03 


1, 


.853 


1, 


.69 


ce.81 






1, 


.934 


1, 


.95 


N6090 


EN2- 


■04 


1, 


.775 


0, 


.26 


ce.91 






1, 


.876 


1, 


.04 


spil.4 


EN2- 


■05 


2, 


.056 


1, 


.17 


122491 






2, 


.077 


1, 


.04 


ce.86 


EN2- 


■06 


1, 


.951 


1, 


.95 


122491 






2, 


.008 


1, 


.95 


ce.91 


EN2- 


■07 


2, 


.152 


1, 


.95 


N6090 






2, 


.339 


1, 


.30 


120551 






2, 


.450 


1, 


.04 


ce.86 


EN2- 


■08 


1, 


.960 


1, 


.95 


N6090 






2, 


.087 


1, 


.17 


122491 






2, 


.174 


1, 


.04 


ce.81 


EN2- 


■09 


2, 


.180 


1, 


.95 


N6090 
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Table 4 — Continued 



Source ID z Av Template 



2.273 1.95 ce_91 
2.299 1.95 120551 



XMM-01 


1, 


.890 


1, 


.95 


ce_81 




1, 


.972 


1, 


.95 


N6090 


XMM-02 


1, 


.928 


0, 


.00 


ce_81 




2, 


.273 


0, 


.78 


122491 


XMM-03 


2, 


.113 


0, 


.00 


N6090 




1, 


.413 


0, 


.00 


N6090 


XMM-04 


2, 


.124 


1, 


.95 


ce_81 




2, 


.212 


1, 


.95 


N6090 


XMM-05 


1, 


.010 


1, 


.82 


SBSey 




1, 


.147 


1, 


.17 


ce_86 




1, 


.778 


0, 


.13 


N6090 


XMM-06 


1, 


.010 


1, 


.95 


M82 




1, 


.205 


1, 


.82 


ce_86 



Note. — The solutions are listed in order 
of . Al l the 'ce' templa t es are from the 
library in ICharv fc Elbaj feOOlh. the oth- 
ers ar e from the library in iPoUetta et af] 
d2007l ). The fits are shown in Figures (TS] 
[141 [Isl and [Tel 
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Table 5 

fir and nir luminosities and dust temperatures for the lockman hole sample 



Source 
ID 


T 

dust 


(Lq) 




u Li 
(L 


\) 

1) ) 


LH-01 


26.25 


12.23 (11.94, 12.46) 


11 


93 


LH-02 


31.76 


12.62 (12.39, 12.79) 


11 


69 


LH-03 


38.79 


12.80 (12.62, 12.93) 


12 


28 


LH-04 


41.36 


12.75 (12.59, 12.87) 


11 


91 


LH-05 


31.29 


12.33 (12.33, 12.33) 


11 


74 


LH-06 


29.39 


12.32 (12.32, 12.32) 


12 


07 


LH-07 


30.72 


12.27 (12.27, 12.27) 


12 


09 


LH-08 


32.53 


12.34 (12.12, 12.51) 


12 


57 


LH-09 


35.09 


12.31 (12.11, 12.47) 


12 


15 


LH-10 


34.62 


12.40 (12.20, 12.56) 


12 


25 


LH-11 


29.20 


12.04 (11.79, 12.23) 


12 


03 


LH-12 


31.39 


12.28 (12.28, 


12.28) 


12 


04 


LH-13 


29.58 


<12.22 (11.97 


12.41) 


12 


58 


LH-14 


32.72 


<12.44 (12.22 


12.60) 


12 


09 


LH-15 


40.03 


<12.63 (12.46 


12.75) 


11 


96 


LH-16 


36.23 


<12.63 (12.44 


12.78) 


12 


13 


LH-17 


31.00 


<12.36 (12.13 


12.55) 


11 


86 


LH-18 


33.95 


<12.53 (12.32 


12.69) 


11 


79 


LH-19 


33.47 


<12.66 (12.45 


12.83) 


12 


06 


LH-20 


31.48 


<12.35 (12.12 


12.53) 


11 


79 


LH-21 


38.42 


<12.64 (12.46 


12.77) 


12 


21 


LH-22 


28.06 


<12.36 (12.10 


12.57) 


11 


98 


LH-23 


22.74 


<11.63 (11.63 


11.63) 


10 


76 


LH-24 


25.31 


<12.06 (11.76 


12.30) 


11 


54 


LH-25 


33.67 


<12.50 (12.29 


12.67) 


12 


58 


LH-26 


31.29 


<12.31 (12.31 


12.31) 


11 


75 


LH-27 


32.53 


<12.57 (12.35 


12.75) 


11 


92 


LH-28 


40.22 


<12.87 (12.70 


12.99) 


11 


93 


LH-29 


33.95 


<12.56 (12.35 


12.72) 


11 


44 



Note. — L^fj^ is the integrated lu minosity between 
42.5 and 122.5Atm llHelou et al.lll988t) obtained using a 
greybody model with T^„g4=3.9+9.5x (l+z), and /3=1.5 
(see text). The luminosities have been obtained at the 
redshift z of each source. In brackets are the luminosities 
derived at z+ = 2+0.12X (l+z) and z- = 2— 0.12x (l+z). 



Table 6 

Results from the co-added MIPS images for the Lockman Hole sample 



Sample 


N 


< /24 > 


/24 


/to 


/l60 






detections mean 




stack mean 










(mjy) 






All 


27 


0.83±0.31 


0.93±0.08 


5.7±1.5 


39.0±10.2 


f 1.2mm > 2(7 


12 


0.92±0.33 


0.99±0.11 


5.7±1.3 


57.8±16.3 


f 1.2mm < 2cr 


15 


0.75±0.28 


0.90±0.11 


7.3±2.8 


16.9±8.2 
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Table 7 

Stellar masses for the Lockman Hole sample 



Source ID 


Av 


(10" Mq) 







■Ul 


i.2ii«: 


07 
03 





oo + O, 

.32_o 


,02 
,03 


LH- 


■02 


i.5il'>: 


14 
,06 


2, 


,27li', 


,05 
,18 


LH- 


■03 


2.53+;;' 


,36 
,18 


1, 


,7011', 


,06 
,10 


LH- 


■04 


4.66+;!' 


19 

,1.'5 


2, 


•Sill]' 


,03 
,09 


LH- 


■05 


2, 


,64«: 


,39 
,03 


1, 


.851- 


,02 
,15 


LH- 


■06 


1, 


,7iir,: 


,08 
,05 


2, 


.371° 


,04 
,06 


LH- 


■07 


1, 


,06i:;: 


,05 
,03 


2, 


.301- 


,03 
,05 


LH- 


■08 


4, 


,94i:;: 


10 

,56 


1, 


,671- 


,02 
,16 


LH- 


■09 


3, 


,25i;;: 


25 

,14 


1, 


,801- 


,04 
,09 


LH- 


■10 


1, 


.6iis;: 


10 
,05 


2, 


,131- 


,05 
,08 


LH- 


■11 


2, 


.3411!: 


15 
,04 


1 


.931- 


,03 
,06 


LH- 


■12 


3, 


,i4i;;: 


31 
,20 


1 


.251- 


,03 
,04 


lh- 


■13 


1, 


,43i:;: 


29 
,12 


1 


.761- 


,02 
,17 


LH- 


■14 


1, 


.82i;;: 


15 
,08 


1 


.851- 


,04 
,16 


LH- 


■15 


5, 


,6oi:;: 


15 

,14 


1 


.861- 


,05 
,05 


LH- 


■16 


1, 


,931- 


17 
,09 


2 


.731- 


,03 
,13 


LH- 


■17 


2, 


,701": 


,31 
,15 


2 


.051- 


,04 
,10 


LH- 


■18 


0.85+0' 


12 
,03 


2 


.851- 


,02 
,20 


LH- 


■19 


2.31+;!' 


23 
,13 


1 


74+0, 


,20 
,22 


LH- 


■20 


1.431;;; 


33 
,09 


1 


.791- 


,16 
,36 


LH- 


■21 


4.371": 


27 
,12 


1 


.881- 


,05 
,06 


LH- 


■22 


0.601": 


04 
,02 


2 


.491- 


,08 
,15 


lh- 


■23 


0.231": 


06 
,00 


2 


.821- 


,02 
,29 


LH- 


■24 


0.90+" 


07 
,03 


2 


40+0, 


,03 
,04 


LH- 


■25 


1.531- 


25 
,04 


2 


.171- 


,02 
,19 


LH- 


■26 


O.81I;;: 


08 
,02 


2 


.761- 


,03 
,24 


LH- 


■27 


2.22+0 


51 
,03 


1 


.421- 


,02 
,21 


LH- 


■28 


1.86+g' 


,27 
,10 


1, 


.761°: 


,05 
.09 


LH- 


■29 


1.491°: 


,22 
,06 


1, 


.981°: 


,16 
.26 



Note. — Uncertainties correspond to 

1(T. 
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